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AN OSCILLATION CRITERION FOR SECOND ORDER 
NONLINEAR DIFFERENTIAL EQUATION 


S. R. GRACE 


Department of Mathematical Sciences, University of Petroleum and Minerals 
Dhahran, Saudi Arabia 


(Received 7 December 1987) 


A new oscillation criterion is given for the second order nonlinear diffe- 
rential equation 


x(t)+ qi f(x()) = 0 


where the coefficient g(t) is not assumed to be nonnegative for all large 
soe 

values of t. Condition on f of the form | a < ec, used by Onose, Philos 

and Wong is discarded. 


1. INTRODUCTION 


Consider the second order nonlinear differential equation 


FO +4dOSEM) =4(.= a) A) 


where q: [to, ce) > R = (— oe, c0), f: R > R are continuous and x f(x) > 0 for 
are A). 


We consider only those solutions of eqn. (1) which exist on [to, co). A solution 
of eqn. (1) is said to be oscillatory if it has arbitrarily large zeros, and otherwise 
it is said to be nonoscillatory. Equation (1) is called oscillatory if all such solutions 


are oscillatory. 


Recently, Philos? considered the strongly superlinear equation of the form (1) 
i.e., the function f in eqn. (1) is such that 


—-c 


+co Pe Fe 
\so- A < e0 and ny <? a2) 


and proved the following oscillation criterion. 


Theorem A—Suppose that condition (2) holds, 


f@) > k > Oforx#0,(’= 3) ...(3) 
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and let P be a positive continuously differentiable function cn the interval [fo, 9.) such 
that P is nonnegative and decreasing on [/0, o- ). Equation (1) is oscillatory if 


i inf | P (s) g(s) ds > —co ...(4) 
ty 
and 
Para! 
ee - | | P (u) q (u) du ds = ©, 1) 


~ 


o 
> 


0 


His criterion extended and improved some of the results due to Onose! and 
Wong® 5. Theorem A is only concerned with oscillatory behaviour of strongly super- 
linear equations of the form of eqn. (1). Therefore, the purpose of this paper is to 
establish a new oscillation criterion for eqn. (1), where condition (2) is discarded. We 
also discuss the asymptotic behaviour of the forced equation 


x (t) + g(t) f(x (1) = e(t) ...(6) 


where e : [to, .) > Ris continuous. 


Our results extend and improve some of the results of Onose!, Philos? and 
Wong? 5. 


2. MAIn RESULTS 


Theorem 1—Let condition (3) hold and Suppose that there isa differentiable 


function 


P :[to, co) = (0, co) 


such that 
f(t) O fort > io and lim sup { P (a) ds <"es (7) 
t ‘ 
r l 
t . 
2 
micrar _ 2 (u) z, 
sen in | [P (s) q (s) 4b O(n} Ps)! ds > —~ co ...(9) 
to 
and 
l a 
Tina enD ar | t (u) q(u)— 2 Pe) ) | du ds = co, ...(10) 
fit 


0 
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Then every solution of eqn. (1) is oscillatory. 


Proor: Let x (t)be a nonoscillatory solution of eqn. (!). Without loss of 
generality, we assume that x (t) #0 for allt > fo. Furthermore, we assume that 
x(t)>0 fort > fo, since the substitution u = — x transforms eqn. (1) into an equ- 
ation of the same form subject to the assumptions of the theorem. 

Now we define 


w (t) = P(t) jo for every t > fo. 


Then for all ¢ > to we obtain 


w(t) = — P(t)q(t) + ao w(t) — XG ——— f’ (x (t)) w2 (ft) 


— 0a —FeO1— Sotvewy — SO p, 





4kp (t) P (t) 2/k 
(11) 
Hence, for all t > ty we have 
t - t 
w () <w(te) — | (009 (9) — SA) as - Je rey Wk w Os) 
Zi) 
ays 7 iP ds. (12) 


Next, we consider the following three cases for the behaviour of < : 


Case |—% is oscillatory. Then there exists a sequence {tm}m.1,2...in [fo, co) 
with lim tm = oo and such that * (tm) = 0(m = 1,2,...). Thus, (12) gives 


tm ° tm 
P (s) p2 (s) 
| a IVE w ©) — SU EP as < w (0) - J (6 (8). 466) — Oo 
t 
7 ariel, 2,, : : 


and hence, by (g) we conclude that 


Ps) ia 
[agrvero - ap ee ae 
t 


0 

So, for some positive K we have 
eo ° 

Jay p(s) [A w(s) — ar }? ds & K for every t > fo. 


ty 
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By the Schwarz inequality, for t > fo we get 


t 


: t 
1— | vk wo fe ya12=1 | ve@l 
0 0 





THEO 


t t 


p(s) es 
TR asl? < (JP@ads | 7 0G) Le Oe 








ty to 


t 
< K f§ P(s)ds. 
t 


From (7), there exists a positive constant Z such that 


t 
J p(s) ds € Lt? for every t > to 
to 


and hence for all t > to we obtain 





t . 
—|tvew@— £Ojas< vert, 
t 


0 


Since P (t) > 0 for t > to we have 


t 
=I w (ds < [RE t. 


Furthermore, (12) gives 
t 


J P96) — Oras < — wi) + w(t 
"9 


and therefore, for all t > to 


t & 


t 
W) 
J fe waw—~ 2 das < < — [| wo) ds + = 0) w(t) 
0 





to 9 t 


Pray | €é 
S Nee t + (t — to) w (to) 


and hence 


~ 


“™ | ee 
ny 
ous 


(u) 
[p (u) q (u) — oom du ds < SE +( 1 — 7) w (00). 


~ 
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This contradicts condition (10). Thus * (¢) is of constant sign for all t > fo. 


Case 2—% > 0 on [h, cc) for some ft; > fo. In this case, from (12) it follows 
that fort > th 


t 4 
te (s) q (sy) — 2) ds < w(t) 





4kp (s) 
to 
and consequently 
t 8 . 
2 
Lf fir eu — F305 aed = Sven 
'y to 


which again contradicts condition (10). 


Case 3—% <0 on[M, cc) for some ti > fo. From (11) and the fact that P is 
nondecreasing on [fo, co) it follows that 


t = t 
ew (few (a) | © (s) 4) - Fes lds + | xa) 
ty ty 
f' (x (s)) w? (s) ds. (13) 


If 


l 
P (s) 





Tt’ (x (s)) w? (s) ds < 00 


pe ty 


then condition (3) ensures that 





¢ ] 
w2 (s) ds < co 
romana 
4 
and hence we can arrive at a contradiction by the procedure of Case I. So, we suppose 
that the above improper integral is infinite. From (13) and condition (9) we derive 
t 
-—w(fe2ct+ | 


ty 


] 
——— w2(s) f’ (x (s)) ds 
P (s) ( 
where C is a constant. Furthermore we choose a fo > 1 so that 


to 


Cc +| ray w2 (s) f' (x (s)) ds = Ci > 0 


ty 
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and then for every ¢ > fo we get 
t 


we ONC + [Zo we OS) O) a]? 


ty 


1 
p (t) 








_ #W)f' &(D) 


4 FO) 
and hence, by integrating over [f2, t], we obtain 


t 





l ae ee F(x (t2)) 
m8 Ci [c +| P (s) WaT: Saas) oe Ie Sf (x(t) © 
* 
Thus 
l ] 
C+ | p (s) w2 (s) f’ (x(s)) ds > Co Fix)’ for all t > fo 


ty 


where Cz = Cj f (x (t2)) > 0. So (13) yields 


*(t)<— Ce fOr every t D> ta 


1 
P(t) 


and consequently we have 





t 
x() < x(n) — C2 5 a ds +> — coast > co 


~*~ 


2 
a contradiction to the fact that x (t) > 0 for t > to. This completes the proof. 


Theorem 2—Let conditions (7) and (8) in Theorem | be replaced by 


P (t) > Oand P(t) < O fort > to, ; ...(14) 
Then the conclusion of Theorem | holds. 
Proor : The proof is similar to that of Theorem | and hence is omitted. 
From the proof of Theorem 1 we see the following results hold : 
Corollary 1—Let the conditions (3), (7) and (8) hold and 


iva) 


p2 (s) 
| 0 Mesa aan (15) 


Equation (1) is oscillatory if condition (4) and (5) hold. 
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Corollary 2—Let conditions (3), (4), (5), (14) and (15) hold, then equation (1) 
is oscillatory. 


Remark 1: Theorem | includes Theorem 3 in Philos? asa special case, and 
since condition (2) is disregarded in our criterion, Corollary | extends and improves 
Theorem A when condition (15) holds. Onthe other hand, Theorem | can be applied 


in some cases in which Theorem A is not applicable. Such a case is described in the 
following example. 


Example 1—Consider the differential equation 


X(t) + (—t-4P2 sin t + 3-3/2 (2 + cost) f(x (t) =0,t > to = > 





2 
...(16) 
where f is any one of the following functions : 
(i) f(x) = cx,x E R force > 0. 
(ii) f (x) = |x| 7%sgnx + kx,x € Rfory > Oandk > 0. 
(iii) f (x) = x log? (e+ |x|),x © Rforp > lI. 
(iv) f &) — xe4l4l, x E.R, Az O. 
We let p (1) = t. Conditions (3), (7) and (8) are easy to verify and 
t . t 
| [p (s) g(s) — ae ds = | a — /ssin s — ae | a 
t wiz 


0 


ll 


Vt(2 + cost) — 2( + }u2 
1 2t 7 

= aK |" cs >vi-2( 5 jn 
1 





i 
= 
3 


where k is as in condition (3), and 





g-8 : t 
] P2 (u) oa = je 
LLL pada) — Fat duas > {vs —2(3-) 
ty to r/2 
l 2s ae! 4. _7 \s;2 
=a” ——] ds = eet 5 ( 5 


7 1) g 7) Ct SY ee ie ieee 
—2( 5} - Y le + Tae oY ). 
Conditions (9) and (10) of Theorem | are satisfied, and hence every solution of eqn. 
(16) is oscillatory. We may note that Theorem A can be applied to eqn. (16) when / is 
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strongly superlinear i. e. for the case (ii) with y > 1, case (iii) and case (iv) with 4 >0. 
Also, Theorem 3 in Philos? is not applicable here since 


8 


fetes sin u 





] du ds < o0, 


2u Vu Vu 


l 
lim sup — 
t-0o se 


nw 


w /2 
We conclude that the class of the function / described is larger than that in Philos?2. 


The following theorem is concerned with the asymptotic behaviour of all solutions 
of the forced equation (6) 


Theorem 3—In addition to the hypotheses of Theorem 1, we let 


fp (s) | e(s) | ds < ©0 method 


then lim inf | x (¢) | = O for all solutions x of eqn. (6). 
t->0 


Proor : Let x be a solution of eqn. (6) on [to, cc) with lim inf | x (t)| > 0. 
—>0o 


Clearly x is nonoscillatory. Without loss of generality, we assume that x (t)> 0 for 
t > to. Furthermore, we consider the function W defined as in the proof of Theorem] 
and then for ¢ > to we get 


PR eas : 
w(t) = p(t) Fou & ~49@ tem) P(t) e(t) — Pit) 


w2 u < = ] a, 
(f(x ())< — P(t) q(t) + t(o) P(t) | e(t) | 


x (t) 


1 
+ p (t) Tey PH w2 (tr) f’ (x (t)) 


where c = inf x(t) > 0. Thus 
t7 to 


P (s) g (s) ds < —w (t) + w (to) + | ES) ww (as 


or a, Ch 


to 
f t 
+ Foe | Pte d l 
Fee JPME ds — | 5 we wy" (xo) ds 
0 ty 


Ce < t 
ae P (s) 
<M— wt) + | FG) ds — loa w2 (s) f# (x(s)) 
ty to ds. 
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co 
] 
where M = w (fo) + ror | P (s) | e(s) | ds. Now, we can complete the proof by 
t 


procedure of the proof of Theorem | and hence we omit the detail. 
The following examples are illustrative. 


Example 2—Consider the second order forced equation 





Aint 2+ cos t sin ¢ 15 ho ee 
x(t) + Mook eae Els (t) = hie fo sin f 
Lor t 
ee eth (18) 


All conditions of Theorem 3 are satisfied ants p(t) = t. Hence, all solutions x of 
equation (18) satisfy ise inf | x(t) | = 0. One such solution is x (t) = I/t/t. We 


may note that peipeation 1 in Philos? is not applicable to equation (18). 


Example 3—The equation 
ees hee td cos t, t > 0 ..(19) 
tg I 
has the oscillatory solution x (t) = an All conditions of Theorem 3 are satisfied 
with P(t) = 1,t >to > 0. 
Remark 2: From Example 2, we see that eqn. (18) without forcing term is 
oscillatory by Theorem | while the forced equation (18) has a nonoscillatory solution 


x(i)= +> 0Oast— co. Onthe other hand, from Example 3, we find that the 


unforced equation x +x =0is oscillatory and has solution sin t and cost * O as 
t > co and eqn. (19) has oscillatory solution x (t) = sin t/t > 0 as t + oe. Therefore, 
it seems that the size of the forcing term is responsible for generating such behaviours. 
Now, it remains an open problem to find conditions on the forcing term “‘e” in eqn. (6) 
so that the oscillatory character of the unforced equation is either to be changed or 


maintained. 
Remark 3: The results of this paper can be extended to more general equations 

of the form 

(a(t) #(t)) + p(t) #(t) + g(t) f (*(t) = 0 ...(20) 
and 

(a(the(t) + pMae(O+ a) F&M) =e .. (21) 
where a, e, p,q: [to, 00) + Randf: R > R are continuous, a (t) > O for t > to and 
x f(x) > Ofor x #0. 
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ON THE UNIFORM STABILITY OF A SYSTEM OF DIFFERENTIAL 
EQUATIONS WITH COMPLEX COEFFICIENTS 


Z. ZAHREDDINE 


Department of Mathematics, U.A.E. University, Al-Ain, P. O. Box 15551 
United Arab Emirates 


(Received 9 June 1988) 


In a recent paper‘, necessary and sufficient conditions for the asymptotic 
stability of a system of differential equations of dimension at most 4 with 
complex coefficients were established. In the present work, other types of 
stability are considered, the results of Zahreddine and Elshehawey* are 
extended, and necessary and sufficient condition for the stability and uniform 
stability of the above systems are established, 


1. INTRODUCTION 


Consider the homogeneous, first order linear system of ordinary differential 
equations of n-dimensions X’ = AX where A isann xX nm complex constant matrix 
and ¥ (t) is a column vector of the n dependent variables. The characteristic 
equation is a polynomial equation of degree n vhose roots are the eigenvalues of A. 


We follow the definitions of asymptotic stability, uniform stability and stability 
as given in Jordan and Smith?, and we note that when 4 is constant, then stability im- 
plies uniform stability [Jordan and Smith?, remarks following definition (9.3)]. 


By (Jordan and Smith?, Theorem 9.3), the question of stability of the system 
X’ = AX is related to the nature of the eigenvalues of A, when A is complex, more 
details may be found in Boyce and Diprima!. With the help of Hurwitz polynomials 
and positive functions, we were able to establish (Section 3 of Zahreddine and 
Elshehawey4) necessary and sufficient conditions for the asymptotic stability of the 
system X’ = AX where A is a complex matrix of dimension at most 4. We intend to 
include not only asymptotic stability, but also stability and therefore uniform stability. 
We extend the results of Zahreddine and Elshehawey? and establish necessary and suffi- 
cient conditions for the stability and uniform stability of the systems described above. 


Before we proceed, we recall some basic definitions and facts. 


2. DeFINITIONS AND NOTATIONS 
Definition 2.1\—The polynomial f (4) = A” + 41 An-1 + ... + an-1A + an with 
complex coefficients is a Hurwitz polynomial if all its roots have negative real 
parts. 
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Definition 2.2—If g (A) is a rational function, its paraconjugate is defined by g* (A) 
= Ete. where A denotes the complex conjugate of A. 
When f (A) =A" + a, A™1 4 ...... + an_2 2 + an_1A + an 
then ; 
S* (A) = (-— 1)" A" + (—1)™1 ay AM FO. + ang A2 — Gn_y A + aan, 
Definition 2.3—A function h (A) is said to be positive if Re A (A) > 0 whenever 
Re A> 0. 


The study of Hurwitz polynomials may be reduced to that of positive functions 
Levinson and Redheffer®, (Theorem 5.1), and, by Theorem 5.2 of Levinson and 
Redheffer? any rational function A such that h and —h* are both positive can be 
written in the form: 

be bn 


by 
h(d)=a+bdA4 \—in, > Aosiwgr et eae sa Ch) 





where Rea = 0,b > 0, bk > O and where the Wj are distinct real numbers. This 
form will be referred to very frequently in the arguments that follow. According to 
the proof of Theorem 5.2 in Levinson and Redheffer?, the W; are the roots of 1/h and 
since they are distinct, it is easy to show that the expansion (1) of A (A) is unique. 

3. STABILITY OF THE SystEM X’ = AY 


We need the following two lemmas. 


Lemma 3.1—If f and f* have only one root in common, then it must have 
zero real part. 


PROOF : Write f(A) and f* (A) in the factored forms : 


f (A) = (A— Aq) (A = Ag)... . (A — An) 
and 


f£* (A) = (-1)" (A+ Ag) (A + Ae) ...... (A + An). 


Suppose that A is the only common root to f and /*. 


If Ay = — dj for some J between 2 and n, then j= — y. Since — ‘i is a root 
to f*, As becomes another common root to fand f* leading to a contradiction. 
Hence Ay = — A1, or Re Ay = 0, 


Lemma 3.2—A1 isa common root 


tofand f* if and only if it j 
root to f+ f* and f — f*, if nd only if itis a common 


Proor : It follows from the obsery 


ion that, f = } * st 
and f* = 31+ 7°) — (f— fey, ae 
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Now, we consider the case where A is a 2 x 2 complex matrix with characteri- 
stic polynomial 


fA) =M+aAt a. Let f*(d) =’2— A+ a 
be the paraconjugate of /. 


Theorem 3.1—The system X¥’ = AX where X isa 2 x 2 complex matrix with 


no repeated zero eigenvalue, is stable if and only if one of the following two conditions 
hold : 


1. Rea, > Oand Rea; Re (a; az) — (Im ag)? > 0 


2. Re ay = Im ag = 0 and aj — 4a < 0, where asymptotic stability occurs only 
when both inequalities in 1, are strict. 


Proor : The case of asymptotic stability is settled by Theorem 3.1 Zahreddine 
and Elshehawey?. 


Uniform stability, or stability which is not asymptotic occurs only in each of the 
following cases : 


Case 1— f (A) = (A — i Ax) (A — Ag) where Re Az < 0 

Case 2— f(A) = (A — i A1) (A —iA2) where Ag isa non-zero real number. In 
both cases A1 is a possibly zero real number? (Theorem 9.3). 

Case 1— f(A) = (A — 13) (A — Ad). 
Define 


We) a a 


Since Re Ao < 0g (A) is a Hurwitz polynomial, 


f* 0) 


g* (A) = — ee = — (A + do). 


Let 


a BA) =s* 0) fA +f" 
aa g (A) + g* (A) fA) —f*Q) 


g and g* have no roots in common! (Lemma 3.1) 


Therefore h is a positive function? (Theorem 4} 


Since h* (A) = — A (A), — h* is positive, and h (A) can be expanded uniquely as in (1), 
(Levinson and Redheffer?, Theorem 5.2) 
2 + Im az A + Re ao, 2A + Ag - Ag _ 2.02) 


aero ott imo. OSC (Ag + Ad 
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If Re a; = 0, then A (A) becomes a seeond degree polynomial in A which obviously 
is not true. Hence Rea; ~ 9. We execute a long division to bring A (A) to the form 


i 1 


Im a2 a 
Re a2 + TRege > (Re a1 Im a1 — Im as) 


a Re a4, A +ilmage 


which when compared to (1) leads to Re a1 > 0. 


Since / (A) is a first degree polynomial in A, the remainder of this division must be 
zero, implying that Re a; Re (a; 42) —- (Im ae)? = 0 


Conversely, assume that Re a; > 0 and Re a; Re (ay a2) — (Im a2)? = 0. Take 
h (A) as in the first of the proof : 


hay= LOT _ M+ ilmay d+ Reag 
fA) —f* A) Re aj A + i Im ao 


i ] 
eS) «3 Im — Im are ee eee 
(Re a})* 0.4) = Baca (Re a1) A -- (3) 
for the remainder of this division is zero. 


This shows that f + f* and f — f* have a common root, thus impiying that 
fand f* have a common root (Lemma 3.2) which must be unique, for otherwise 
f = f* leading to Re aj = 0. Therefore the common root tof and f* has zero real 
part (Lemma 3.1) call it i Ay, where A; real. 


Let 
f (A) = (A — 7x) (A — Ag) 
since Re a; > 0, by (3), A (A) is positive function. (Theorem 5.2). 
Let 
qe AU 


r = id = A — Ao 
then 
* Sf * (A) = 
g*(A)= - SRST TS = — (A + Ag) 
and 


h(a) = —& A) — 8* (a) 
g(A) + g* (a): 


Now g and g* have no common roots, a 


: nd A positive. : ; 
nomial (Levinson and Redheffer?, FT eee ree ta Plc lie pols 


Theorem ora S leading to Re jo < 0 
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Case 2— f (A) = (A — idx) (A — ig). 
f(A) = Az — i (Ay + Ag) A — Ad Ag = 02 + LAF a2 
so Ay + Ao = i aj and Aj Ag = — azo. 


Hence Re a; = Im aa = 0. it is easy to check the identity (A1 + A2)? — 441 A2 = (Ai 
— A? > 0 


Therefore 4a2 — aj > 0 
conversely, assume Re a; = Im ag = 0, and 4az — aj => 0. 


Consider the quadratic equation with real coefficients X 2 — ja, X¥ — ag = 0. Its dis- 

criminant — a? + 4a > 0, so it has two real roots A1 and A> such that: 

Ay + Ag = ia, and Ay A2 = — ag. It is easy to verify that f(Q) = (A — id1) (Aida). 
Consider now, the system ¥ = AX where A isa3 X 3 complex matrix with 

characteristic polynomial 


f A) = A3 + a1 2 + aA + a3, and let f * (A) = — 8 + a A2 — a2 A + ag be the 
paraconjugate of /. 


Define the numbers a, b and c by: 
a — Re aj Re Re (ai 42 — a3) — (Im az)’, 
b = Re ay Re (ae a3) — (Re a3)”, and 
c = Rea, Im (41 as) + 43 Im ap. 


Theorem 3.2—The system X’ = AX where Ais a 3 x 3 complex matrix no re- 
peated zero eigenvalue, is stable if and only if one of the following three conditions 


hold : 

if Rea, > 0, a > Oand ab —c? 2 0 

2. Rea, > 0, a=c = Oand 4 Rea Re (a1 42) — 3 (Im am)’ > 0 

3. Rea, = Imag = Re a3 = 0 and there exists a number « with Re « = 0, f («) 


= 0 and satisfying 3 «? and 2 a1 « — a? + 4ag > 0, where asymptotic stability occurs 
only when all the inequalities in 1, are strict. 


Proor : The case of asymptotic stability is settled by (Zahreddine and Elshehawey?, 
Theorem 3.2). Uniform stability, or stability which is not asymptotic, occurs only in each 
of the following cases: 


Case 1—f (A) = (A—id1) (A - do) (A — As) where Re Ao < 0, Re A3 < 9. 
Case2— f(a) = (A — iA1) (A — ide) (A — As) where Ao is a non zero real 
number and Re Az < 0. 
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Case 3— f (A) = (A — idx) (A — ida) (A — idg) where Ag and Ag are non-zero 
real numbers. 


In all three cases, Ai denotes a possibly zero real number (Theorem 9.3 of 
Jordan and Smith2). 


Case 1— f(A) = (A — idx) (A — Aa) (A — Az) 
Define 


gQ) = LO = @ -aya— rs) 


g is a Hurwitz polynomial for Re A2 < O and Re A3 < 0 


* = POE = AW? i 
g* (A) = — pene = (A + Aa) (A + Ag) 


let 


ha) = 2M—2°@ _ SWM+f*O) 
8 (A) + g* (A) f(A) —f* @) 


g and g* have no roots in common (Lemma 3.1 of Zahreddine and Elshehawey4) and 
therefore / is a positive function, (Theorem 5. | of Levinson and Redheffer?). 


Easy to see that 


h(a) = Re a1 A2 + i Im aed + Reas 
i AS +i Im a; A2 + Re aod + i Im a3 


For any complex number y + 0 it is easily checked that Re y and Re (I/y) have the 


Same sign. Hence 1/h is positive if h is positive. And since h* (A) = — h (A), then 
(1/h)* = — I/h, hence — (1/A)* is positive. Therefore the function 

| — A®° + iimay A2 + Reaga + ilImas 4 

h (A) Re aj A? + i Imas A + Re ag (4) 


can be written as in (1), (Theorem 5.2 of Levinson and Redheffer3) 
Also, 
I A* — i Im (Ag + Az) A+ Re (Ag dg) 
h(a) ~ “= Re (to + As) A + i Im (Ag A3) ‘a7 
where Re (Az + Ag) + 0, 


By comparing (4) and (5) we conclude that Re a1 ~ 0. In (4) we execute a long 
division : 
I 


ep heen l 
ROT (Rea? (Re a; Im a, — Im a2) + 


Re ay 





A 


ob Z 
Do eee 
Re ay 2 + i Im az X + Re ag ‘ 
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where 


l 
Ae ae! Pe \ 
R ea (aA + ic) 


a and c are defined above. 


Re a; > 0 (Theorem 5.2 of Levinson and Redheffer3). Also Theorem 5.2 of Levinson 
and Redheffer? and the remarks that follow as illustrated by Example 5.1 in Levinson 


; A2 [ ‘ : $4 
and Redheffer? imply that K (A) = peat eee tn Gade -t Re dg is a positive 


function, and so is — K* (A). Therefore K (A) has form (1), where again it can be 
shown that a/(Re a,)? the coefficient of \ in R is non-zero. 


By (5) 1/h reduces to a second degree polynomial in A over a first degree one, 
which in turn implies that 


A2 
Rea 0 + imag + Re as appears as a first degree polynomial in A and 
(2) of case 1, Theorem 3.1 has re-emerged but now with different coefficients. A 
repetition of that argument leads to a > 0 and ab — c2 = 0, 


Conversely, assume that Re a, > 0,a > Oand ab — c2 = 0 


f(A) — f* A) se A3 + i Ima, a2 + ReagrA +ilmag (6) 
J (A) + f* (A) Re a1 A2 + ilmaga + Reag nee 
i l 
= (Re ai)2 (Re a; Im a; — Im a2) + Re a1 A 
R 
+ 


Re a; A2 + iImagiA + Reags A 
where R as defined before. 


R 


= ha . : SN a eee 
ab — c2 = O implies that Re a, A2 + ilmaz\A + Reag 


reduces to a non-zero constant over a first degree polynomial in A. Therefore f + f* 
and f — f* have acommon root which must be unique, for otherwise (f — f*)/( f 
+ f*) reduces to a first degree polynomial in A, leading to R = 0 which is certainly 
not true since a > 0. Therefore f and f* have only one root in common (Lemma 3.2), 
which must have zero real part (Lemma 3.1), call it 7 A1. 





Let 
SQ) = A — idx) (A — Aa) (A — Asa). 
Now if 
(A) +f*@) 1 | fM=f*@ 
oA) a ae — f* (A) fee aay fa) +f*A) 


is a positive function since Re a1 > O anda > 0. So his also a positive function. 
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Let 

£0) = LO -@- 9 a-29) 
then 

ery = - AY - 0 + Wat WD) 
and 

Kays g (A) — g* (A) 


g (A) + g* (A) 


g and g* do not vanish simultaneously since i A is the only common root of f and f *, 
and because A (A) is positive we conclude that g (A) is a Hurwitz polynomial (Theorem 
5.1, Levinson and Redheffer®) therefore Re 42 < 0 and Re A3 < 0 


Case 2— f(A) = (A — iA1) (A— ida) (A — As) 


oo D(A) So eo 
.. Oa Ue 


is a Hurwitz polynomial since Re A<0 


g (A) = (A — tax) (A — 7 Ag) = —(A+ Ag). 


Let 
h OV me = Se LOE FeO) 


g(y+ae*() f(A) + f* (a) 


h is a positive function (Lemma 3.1 Zahreddine and Elshehawey* and Theorem 5.1 of 


Levinson and Redheffer*), and sois — h* = h therefore h(A) takes from 1), (Theorem 
5.2 of Levinson and Redheffer?) 


Also A can be written : 


h(a) = + ilmay A2 + Reazr+itmas _ 2 — (As — 2s) 
Rea, 42 +ilImacsA + Rea x = . 
. ~(Ag + As) 


| tea 
It is easy to see that Re a; + 0. The first part of (7) leads to 


i 
h (A) = (Re a1)2 (Re ay Im ay = Im ae) + 





Re ay A 


R 
SS 
Re a, 4 + iJm ao A + Re ag 


x i l : 
By (7), R=0. Since R = (Re a2 (a\ + ic), thna=c = 0. 
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Re a; > 0 for Ais positive (Theorem 5.2 of Levinson and Redheffer®). 
[t remains to verify the inequality 4 Re a; Re (a1 42) — 3 (Im ag)? > 0. 
Since 
f (A) = (A — i Ax) (A — 7 AQ) (A — Ag) = A3 + ay 2 + a2 A + a3 
then 
iAy tide + Az = — ay, iA1 Ag + i Ag AQ — Aq Ag = a, A1 Ag AZ = Az. 


The real part of the first relation, imaginary parts of the second and third lead 
respectively to : 





Re A3 = — Rea, (Ai + Az) Re A3 = Im ag and A Ag Im Ag = Im az. 
Therefore 
Im ae Re az 
‘ = an > Saar ame A ae 
My + Ag Rea and Aq A: Ree, 


So A, and Ag are the roots of the quadratic equation 


Re a; X¥2 + Im ap X — Re az = 0. 


Since Ay and Az are real (Im a2)? + 4 Re ai Re ag > 0, but a = Re a; Re (a1 ag — a3) 
— (Im a2)? = 0 which implies (Im az)2 + 4 Rea; Re a3 = 4 Re ay Re (a1 42) — 3 
(Im az)?, leading therefore to the conclusion. 


Conversely, suppose that Re aj > 0,a=c= 0 and 4 Re ay Re (a; 42) — 3 
(Im ag)? > 0 


“OEY Ech ile €,9 immer ed + iIma,A?2 + ReagA +1imag (8) 
Fi) + ft A)” Re a; A2 + i Imag A + Reag 2 


(Re a; — Im a,;—Imag) + 


ae (Re ay)? Re a1 


for the remainder of this division is zero, sincea = c = 0. Therefore f — f* and 
f +f * have two common roots, implying that f and f* also have two common roots 
(Lemma 3.2), call them Ai and 9. 


Let 
f A =(Q— AQ = 42) A — As) 
As cannot be another common root to f and f*, for then f= —f* leading to 
Re ai = 9 
write 


ft) =—-AtAW@A + Aa) (A + 3) 
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if 41 = — As then Az = — A, ands ince — Aj is a root Of f*, Az becomes a common 
root to f and f * which is impossible. The same applies on Ag. 


Therefore we have only two possibilities either A1 = — AL leading to 
Aa, = — As 
orA,; = — de which is equivalent to Az = — oe 
In both cases, we get f* (A) = — (A — Ax) (A — Ag) (A + As). 
Let 
f(A) 
ReShaciie  egieea a) 3 
so 
‘je A Oe ee 

GPA Sa eet : 

Let 


& (A) + g* (A) FA) ASE) 


In (8) above, Rea; > Oand A* = —h therefore h (A) is a positive function 
(Theorem 5.2 of Levinson and Redheffer®), and since g and g* have no roots in com- 


mon, g is a Hurwitz polynomial, (Theorem 5.1 of Levinson and Redheffer?). Hence 
Re A3 <0. 


hoy = £M=2@ _ ff 


—_ { 
[t remains to determine the nature of \; and Ag. When Ay = — Ay andAg = — dos 
then Re Ay = Reds = 0, 


In this case the inequality 4 Re a; Re (a1 ag) — 3 (Im a2)? 


= 0, is easily verified as in 
the previous section. 


Consider now the second possibility where Ay = — oy 


or Az = — dy. This leads to Re 41 = — Re Av and 
Im 4; = Im Ag. Call B = Re 41 and y = Im )4, so 
M1 = 8B + iy, andre = — 8 + iy, 
fM=-PA—(@+ HP -(—p+ iy)]. (A — As) 
= AS + ay 02 + aed + az, Then : 


iy + As = — ay, Diy Ap — B2 — y2 = ay (pe 4 y?) As = ag. 
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Obtain 43 from the first relation and substitute it in the second to get: 3,2 
— 2iayy — B® — ap = 0 whose real and imaginary parts produce : 


27 Rea, + Imag = 0 and B2 = 3y2 + 2 Ima, y — Re az. 
In the last relation, we substitute y = — Im ao/2Re aj, to finally have 
a2 = 3 (Im ag)? — 4Re ay Re (a1 a2) 
4 (Re a)? 


But this requires, since #2 > 0, that 3 (Imag)? — 4Re a; Re (ai 42) > 0, and by 
assumption, we have 3 (Im a2)? — 4Re a; Re (a; 42) & 0, when we compromise these 
two inequalities, we end up with 8 = 0 leading to Re A1 = 0 and Re Azo = 0. 


Case 3— f (A) = (A — idx) (A — ida) (A — ids) 
= A3 + ayd?2 + aad + az. 
Therefore 
Ay + Ag + Ag = iay, Ay Ag + Ag Ag + AZAL = — ag and Aji Ag As 
= — ji a3. 
These relations imply Re a; = Im a2 = Re ag = 0. 


The first and second relations lead to Ag + A3 = iay —A1 and Ag A3 = — a2 - AL 
(Ag + A3) = — ag + A? — ia, A1. 


So Ao and Az are the roots of the quadratic equation with real coefficients X? + (A1 
— jay) X + (A? — iay Ay — de) = 0. 


Since Ao and Ag are real, (Ay — ia1)2 — 4 (Aj — iai Ay — a2) > 0, 
or — 3A2 + 2/a, \1 — a? + 4 a2 > O. Letting « = A1, we get 
302 + 2aja — af + 4aq > 0. 


Conversely, suppose Re a; = Im ao = Re az = 0, and there exists « with Re « = 0, 
f («) = 0 and 3x2 + 2a, a — a? + 4aq > 0, form the quadratic equation with real 
coefficients. 


X2 —i(a +a) X¥ — (a2 + 14 + ag) = 0. The discriminant 32 + 2aja — aj 
+ uao > 0, therefore there exists two real numbers A and A3 such that : 


Ao + Ag = i (a + a1) and Ag. As = — (a2 + ayx + ag) 
also, 

Ag. Ag = — ag —a(a +41) = — ag + ia (Ag + Ag). 
Therefore, 


— ia + Ko + A3 = id and — ixAg ~_ ixAg + Ag Ag = — 42. 
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Since f(«) = 0, then #3 + aja? + ag% + a3 = 0, and « (a2 + aja + az) = — a3 
which leads to — «% Ag Az = — 43 

or & Ao As = ag. It is Dow straightforward to verify that 

f(A) = (A — @) (A — idg) (A — ida) and the proof is complete. 


Finally, we consider the system X¥’ = AX where A isa 4 x 4 complex matrix 


with characteristic polynomial, 
f (A) = A4 + ay M3 + ag A? + ag A+ a4 
and let 
f* (A) = At — a AB + ag A2 — ag A + ag be the paraconjugate of f. 


Let a be as defined before Theorem 3.2, and define the numbers r,s and? as 
follows : 


r =a. [Re ay Re (az 43 — ay aq) — (Re a3)?] — [Im a, Re (a, — a3 a4) 


+ Reag Im ag)}?. 


s = a. [Re ay Re (a3 a4) — (Im ay)?] — [Re a; Re (a1 a4) 


— Im a2 Im ag}. 
t = a. [Rea Im (ae aq) — Re ag Im a4] 


+ [Re ai Re (a1a4) — Im ag Imag]. [Re a, Im (aya3 —a4) + Re ag 
Imag]. 


Theorem 3.3—The system X' = AX where A is a4 xX 4 complex matrix with no 
repeated zero eigenvalue, is stable if and only if one of the following four conditions 
hold : 


(1) Rea >0, a>0, r>Oandrs — 12>0° 


(2) Reay >0,a>0,r=1t= 0 and the inequality 


4 Re a) Re (ai a, ) ~ 3 (Ima; )? > 0 wher a = y + ay 


a, = y> + a1 y + ap, holds whenever Re y < 0 and f(y) = 0. 


(3) Reay > 0,4 =r =0, Rea Re (a a 
x 144) — Imao Ima = 0, a 
exists a real Ay with f (iA 0) = 0 and Satisfying 7 me 





343 +2 Im ae Ay — (ima)? 4Raz 


Re ay Ro a)? Reais 
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(4) Rea; =Im a2 = Re a3 = Im a4 = O and there exist two numbers « and f with 
Re « = ReB = 0, / («) = f (8) = Oand satisfying 3 «2 + 2a,’ « — a, + 4a, > 0 


where 


= B + a; anda, = @+a,8 + ae 


where asymptotic stability occurs only when all the inequalities in (1) are strict. 


Proor : These case of asymptotic stability is settled by Theorem 3.3 of Zahreddine 
and Elshehawey4. Non-asymptotic stability occurs only in each of the following cases. 


Case 1— f (A) 


I 


(A — iA) (A — Az) (A — Az) (A — Aa) where Re Ai < 0, 


pre 93, 4; 


\ 


Case 2—f (A) = (A —iA1) (A — 1 Ag) A — Ag) (A — Aq) where Az is a non-zero 
real number, Re A3 < 0 and Re 44 < 0. 


Case 3—f (A) = (A — iy) (A — ide) (A — iA3) (A — \4), where Ag and Ag are 
non-zero real numbers, Re 44 < 0 


Case 4—f (A) = (A — id) (A — ide) (A — ida) A - iAy), where Az, A3 and Aq are 
non-zero real numbers. 


In all these cases, A, denotes a possibly zero real number (Theorem 9.3, Jordan 
and Smith?). 


Case 1—f (A) = (A — iA1) (A — Az) (A — As) A — Aa) 


g* (A) = ae we ) = (A — Ag) (A — As) A — A4) is a Hurwitz poly- 
nomial since Re At < 0, i = 2, 3, 4. 
g(a) = — ae. — (A + %) @ + Ag) AF Ad). 


Let 


gayt+e*A fA) mat (A): 


g and g* have no common roots (Lemma 3.1 of Zahreddine and Elshehawey*) 
and consequently h is a positive function (Theorem 5.1 of Levinson and Redheffer®). 


M+ itm ay AP + Re ag A2 + iIm ag A + Re aa 
Re aj A® + iim ag A? + Reaza + ilm aq 








h(a) = 


Since h is positive, it takes form (1). Therefore Re a + 0. By executing along diiv- 
sion, we bring A (A) to the form : 
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-A+ ; (Im a, — thy, 
h(a) = “Rea Re ay ‘ Re a1 
R’ 
f 


Re ay A? + +ilm ae dA2 + ReagA+ilmay, 


where 











Re ag Im ae z= Im ag |» 
R = [ Rea: - Ram )+ Rear (Im ay Raat 





R Im a2 
+ if (1m ay - 1m 44 ) _ rer (Im a; — )Ja 


Re a Rea Re ay 
Im aq . Im a2 ,, 
+ Re a4 + Re ay mies Re a1 ) 


Since A is a positive function, Re a; > 0 
Theorem 5.2 in Levinson and Redheffer? and the remarks that follow imply that 


Re a; A® + ilIm ae eee Re agA + ilImay Tae ee 





Now, 





_ g& (A) — g* (A) 
(WSF eats 
= A= As) — As) A — Aa) + (A + Aa) (A+ AB) A+ Ad) 
(A — Aa) (A — As) (A — Ag) — (A + Az) (A + As) (A + Aa) 


It is easy to see that g + g* isa second degree polynomial in A, since the co- 
efficient of A? ing + g* is — 2Re (Ao + As + A4) which is non-zero. 





3 ; 2 ; 
Therefore the positive Acar H+ itm on Af + Re as A+ itm ag should also re- 


duce to a second degree polynomial in A over a degree one. So (4) and (5) of case 1., 
theorem 3.2 reappear but now with different and more complicated coefficients. The 
process here is rather lengthy, but the argument is entirel 


y similar and when repeated 
leadstoa > 0,r > Oandrs — 12 = 0. 


Conversely, assume that Re aj > 0, a> 0, > O-and ret— facet 
Take 


Pioneer A‘ + i Im ay AS + Re az 2 + iImag3aA + Re ay 
f-f* Rear AS +i lmas a2 + ReasA+ilmar 
l 


I 
Re ay A + Rea, (1m an 





repel 


R’ 
eee 
3 Re a1 A? + i Im ag A2 + Reag A +i Im a4 
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where R’ as defined in the first section. If we consider 
3 : 2 ; 
a ea At timo with the set of conditions a > 0, r >0 


andrs — t? = 0, we come to a position entirely similar to (6) in the converse of case 
1, Theorem 3.2, with the same set of conditions but different and cumbersome coeffi- 
cients. Another application of that argument leads to the fact that 


Re a; A? + i Im az A2 + ReagdA + i Imag and R' have only one root in common. 
The same applies on f + f * and /f — f*. By Lemma 3.2 f and f * have only one com: 
mon root which must have zero real part (Lemma 3.1) call it i A1. 


Let f (A) = (A — i At) (A — Aa) (A — Az) (A — Aa). 


If h (A) = j s es F bie , then A is positive? (Theorem 5.2). 


Let 
f(A) sh 
g(A) = Sey kes A2) (A — As) (A — Aa), 
Soh aN = — (A+ Aa) (A + As) (A + Aa) 
so 
pa) = LMS A) _ 2) 8°) 


fA-F*) ge A)+a*” 


Since i Ay is the only common root to fandf*,theng and g* donot vanish 
simultaneously. Hence g is a Hurwitz polynomial? (Theorem 5.1). So Re aAt < 0 for 
i= 2, 3, 4. 


Case 2—f (A) = (A ~ i A1) (A — i Ag) (A — As) (A — Aa) 


f(A) 


g (A) = seta) AA = (A — As) (A — Aa) 


is a Hurwitz polynomial for Re A« < 0:7 = 3, 4. 


g* (A) “es ata Tey = (A oe As) (A + Aa). 


g(a) —s*) _ fa —f*@ 
ga +e* A SATS EA 


h is positive function [Lemma 3.14 and Theorem 5.13] and so is 


Consider h (A) = 


] M4 + i lm a, A? + Re ag A2 + iImagaA_+ Rea 
h ~ Rea, A® + iim az d2 + ReagaA+ ilmas 


1/h takes from (1), and that leads to Re a1 + 


n0F Z. ZAHREDDINE 


Im a2 


Re ay 


] ] 





i 
Py eee ——_——.(| _ 
h Re ay ASH Re ay ans 


Peds 
t Re a, A®° + ilmagA2 + ReagzA+ilmay, 
where R’ is already defined. By (Theorem 5.2 of Levinson and Redheffer®), Re ai > O. 
Also, 


g(a) + e*(A) _ (A— As) (A— Aa) + (A + Az) A — Ad) 


] 
Bo BADEN On) 


The coefficient of A in g — g* is — 2 Re (Ag + Aq) which is non-zero. Hence 1/h 
reduces to a second degree folynomial in A over a first degree one which implies that 
Re a; A? + i Im ag A2 + Re ag A + i Ima, must be divisible by R’. 


We are now in position (7) of case 2 of Theorem 3.2 but with different coefficients. 
The same argument leads to r = t = O, and also toa > 0, since 


Re a; A® + i Im ag A2 + ReagA+ilmay 
Vek 


is a positive (Theorem 5.2 of Levinson and Redheffer3). Let 


Ka) = 8 =a—ava-ada—a) 


= A3 + a, A? + a, A + a, 
where 


4 4 2 ‘ 
a‘ =A4 + 41,4), = rj + a1 As + a2, ay = AY +4123 + az Aq + ag 


But K (A) = (A — idy) (A — ide) (A — Aza) is precisely case 2, Theorem 3.2, and there- 
fore the inequality 4 Re a, Re (a, 4, )— 3 (Im a’, )? > 0 is satisfied. The same 
can be proved with As replacing Aq. 


Conversely suppose Re a; > 0, a > 0,r = t = O and assume the above inequ- 


ality holds with a=y+ 41,0, = yA + ary + ao, whenever, Re y < Oand f (y) 


=e 
Consider 
fesaien _ At +ilmay A3 + Re a2 A? + iImagaA + Reay 
J-—f* Re a, A? + i Im az 02 + ReagA+ ilmag 


(equation continued on p. 323) 
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i oe i Im a2 
ps Re ay Ar Re a, 





R’ 
+ Re aj A? + i imag A2 + ReagA+ilmay, * 

+-(9) 
Now 


acters 2 
Ae A itm dy A+ Rea A + itm aa with r = t = 0 puts 


us back in position (8) in the converse of case 2 of Theorem 3.2 with the same set of 
conditions. A similar argument implies that 


Rea, A? +ilmag A2 + Reg3A+ilmay 


R reduces to a first degree 


polynomial in A. Therefore / + f* and f — f* have two roots in common. By Lemma 
3.2 fand f* have two common roots, call them A, and dg. 


Let f (A) = (A — Ai) (A — Ag) (A — Az) (A — Aa). Ax and Ag are the only common roots 
to f and f*, for if, for instance, A3 is another common root to f and f *, sof + f* and 
* 


f — {* would have three roots in common, which implies shot ee per becomes a 
first degree polynomial in which is not the case. 

Consider f * (A) = (A + Ax) (A + Aa) (A + Az) (A + Ad). 

If Ar = — Ag, then Az = — Aj, and since — dj is a root of f* (A), As becomes a com- 


mon root tof and f* which we proved to be impossible. Similarly Ai cannot 


equal — Aa. The same applies on Ag. Therefore we have only two possibilities : either 


a4 =- A leading to Ag = — Ao, or A1 = — Ag which is equivalent to Ag = — 1, 
In both cases, we get f* (A) = (A — Ai) (A — Az) (A + As) (A + Aa) 
f (Aa) 
= — - =(A— As) (A-A 
Let g (A) eae (Aa Agy.f 4) 

so 

(A + Ax) (A + Aa) ; 
Let 


SAS — 2 (Al f(a) —- f* A 7 1 
g (A) + a* (A) fA)+f* A)’ h (A) 


fA) +f* (a) 
SY - f*W 





h(a) = 
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In (9) above, since Re a; > 0, a> O and (1fh)* = — 1h, 1fh is a positive 
function (Theorem 5.2 of Levinson and Reddheffer?) and because g and g* have no 
roots in common, g is a Hurwitz polynomial? (Theorem 5.1). Hence Re A3 < 0 and 
and Re Aq < 0. It remains to determine the nature of Ai and dag. 


When Ay = — Ar andAe = — Acthen Re Ai = Az = O. 


In this case, the given inequality is easily verified as in the previous section with both 
Ag and Ag. 


Consider the second possibility where A; = — ‘Ag or Ag = - AL. 
Here we have Re Ai; = — Re Ag and Im Ai = Im Ag. 


since Re Ay < 0 and / (Aa) = 0, the inequality, 


4Re a’, Re(a‘, a, ) — 3 (Ima‘, )? > 0 where a) = Aa + a1 


a’, = A, + aida +az holds. 


Let 
(A 
KA) ey a en Ce 
A Aa 
= AP ay At Pd, A+, 
where 
a, = Ay + a AL + avrg t a3 = — a 


If we take up the argument towards the end of the converse of case 2 of Theorem 3.2, 
we can show that Re Ay = Re Ay = 0. 


Case 3— f (A) = (A — i Ai) (A — i Ag) (A — i A3) (A — Aa) 
Fl f * (A) 


(A ar) (ee ae Le A — Aais a Hurwitz poly- 
nomial for Re Aq < 0. 


. j= See AD ree = = 
ew (A — i Ax) (A — i Ag) (A — i As)) = (A fAg). 
The function h (A) = —% (A) — g* (A) 


g(a) + g* (a) is positive (Lemma 3.14 and Theorem 5.13) 
h (A) can also be put in the following forms : 


Ais At + i Im ai 18+ Rea A? + iIma3A + Re ag 
Re ay AS + i Imag A2 + ReazsA+ilmag 


(equation continued on p. 325) 
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_ 2+ raM—As 
— (Aa + Aa) ° 
Since h is positive, then (1) implies Re ai + 0. 
Because h/ appears as a first degree polynomial in A, the remainder R’ of the division 
of A4 + i Im ay A8 + Re ag A2 + iImagA + Re ag by Re a; A2 + i Im ag A? + Re 


az A+ i Im ag mast be zero. That simply implies a = r = U. and Re ay Re (a1 a) om 
Im ao Im aq = 0. 


But 


1 i Im a2 
N =) ee ees sees Eas & | ey) ean 
h <A) Re aj * oie a Rea 


shows that Re a1 > 0, for A is a positive function. 
Define 


Ka) = LO = a —iama—iaya- rs) 


= 8 + a’, 2+ a, A+ ag 
where 


‘ ‘ 2 ‘ 3 2 
Gy = Aa 01,4, Ay + a1 Aa + 42, 0, = Ay T ay Ay + a2 Aa + 93 


Il 


a4 
v4 


= eo 


with K (A), we are in the position of case 3, Theorem 3.2. Therefore the following 


must be true : Re a, = Im a, Re a, = 0, and there exists a real Ag satisfying 


K (i Ao) = 0, hence f (i Ao) = 0 and 325 — 2ia, A+ a‘ — 44, < 0, this inequ- 
ality is obtained by letting « = i Ao in 3, Theorem 3.2. 


Now Rea, = Re (Aa+a1) = 9, implies, Re Aa = — Re 41, and 


Im a, = 2 Re Ag ImAg + Rea Im 44 + Im a; Re Aa + Im az = 9 


Im a2 I 
when we substitute Re Aa = — Re ai we get Im Ag = 7 ee m a4. 
Therefore 


_ Im ae 
= j——. 
1 Re ai 
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Also by substituting the values of Re Aq and Im Aq in a, 


we get 


2 eae ~ Rea ye [Re a1 Re (ai a2) — (Im ap)?] 


and since 


= Rea (Re ai — a2 — az) — (Im az)? = 0 a, = Re asz/Re aj. 


Therefore Ao satisfies 


_Im a2 ae: (Im ag)? Re ag 0 


“Rea, “° (Re ai)? Rea 





3A8 5 
Conversely, suppose Re aj > 0,a = r = O and 


Re a1 Re (a; a4) — Im a2 Im ay = 0 and assume the existence of real Ao with F (ido) 


ee Im ae (Im ag)? Re ag 
= A eee —_ ———_  — 0. 
O satisfying 349 +2 a, (Re a1) Ee 





Take 


f+ f* _ At +ilma A’ + Rea, dA2 +i Imag a + Reag 
f—-— f* Re a; A® + i Imag A? + ReagaA + ilm ay 





the remainder of this division is zero since a = r = 0 and Re ai Re (a; a4) — Im ag 
Im aq = 0, as is easily checked. 


Therefore f + f* and f — f * have three roots in common. By Lemma 3.2 f 
and f * have three common roots Aj, Ag and Az. 


Let 
f(A) = (A — Ax) (A — Ag) (A — Az) (A -— Ag) 
and 
f* A= A+ AA + As) (A + As) (A — Ay). 


It is obvious that Aq cannot be another common root to f and f*, for then f = tee 
leading to Re a; = 0. 


Since by assumption i Ap where Ao real, is a root of J (A), then it must be a com- 
mon root to f andf*. Therefore i Ao is different from Aa. Suppose A1 = 7 Ap. 


f (A) = (A — ido) (A = Aa) (A — Az) (A — Ay) 
and 


J* (A) = (A — i Ao) (A + Aa) (A + As) (A + Ad). 
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If Ao = — Aq then Ay = — Ag and since — Az isa root of f* (A), then Aq becomes a 
common root to f and f* which is impossible. So Ap # — Aa and similarly A3 + a5 Fe 


It is easy to show that we have only two possibilities : 


either Ag = — re leading to 43 = — As, or Az = —A3 which is equivalent to As 
= Ae. 


In both cases f* (A) = (A — ido) (A — Ag (A — As) (A + Aa). 


Let 
Mu f(A ig, es 
oa aa Ooo) 
then 
tee: f* A) one 
SN gata) (A = Aa) (A) eda 
Define 


orp cegee VACA iis lV Peage aC Use aH 
ety ye (0. £ Ait g(a) 


Since Re a; > 0, 


Im a2 
(Im ai ne as ) 





é l i 
b= Re ay Av Re aj 


is a positive function® (Theorem 5.2), and since g and g* have no roots in common, & 
is a Hurwitz polynomial? (Theorem 5.1) implying that Re Aa < 0. 








Also, 
— g(ay— eta) -2A+Ag— As _ gj tas 
he ees g (A) + g* (A) ~ =(Ag + 4a) Reva Re A4 
Hence, 
Im A4 l _ Ima 
Re As = — Re ai and Mey Ree (Im a1 Roe ) 
leading to 
Im a2 
Bas eth | } 
Im A4 oon m a} 
Let 


KA) = LM 2a -iaa-wA-r 


=\3 +a, 2 + a, A+ 4, 


328 Z. ZAHREDDINE 


where 


: Matai, a, =X + aA tas, a, =A, + AS + ao Aa + a3 


2 
i] 


ee ' 
A4 
With the above values of Re Ag and Im Ag, and with the use of relation Re a; Re 


(a1 a4) — Imag Im ay = 0, it is easy to show that Re a’, = Im a, = Re a:, = QO, 
Let « = i Ao, then Rea« = 0 and K («) = 0, Now the given inequality takes the form : 


Im a2 _ (Im az)? Re a3_ 


eo Pate! Rea; (Re a1)? Re ay “a 


But we know from the previous section that : 


a, = i Im a2/Re ai, and a, = Re a3/Re a, therefore 


302+ 2aix— a’ + 4a, > 0. By 3 of Theorem 3.2,Re Az = Re As 
= 0. 


Case 4—f (A) = (A — iA) (A — i Ag) (A — iAzs) (A — ida), Ax is a possible zero 
real number, 2, 3 and 4 are non-zero real numbers. 


Let 


(0) = FO 2 a — inya— ia) Q = ry) 


me Oy AA rt GA + a where 


ree 2 ; 
av=idet+a,a,=—A, + ia, Ag + ag and 
<9 , a 1 
a= — 15 —aA +iad+a,3=- “= <4 
2 2 ; 
i Ag A2 


with g (A), we are in the position of case 3 of Theorem 3.2, therefore Re ay = Im a? 


= Rea, = 0 which simply lead to 


Re ai = Im a2 = Re ag = [may = 0. 
Theorem 3.2 also implies the existence of a number «, with Re « = 0, g («) = 0 
hence f («) = 0 and satisfying 342 + 2 Ua ie nin 4a’, > 0, obviously « coincides 


with either i Ai, i Ag or 7 Aq. If we let B= Az, then Re B=0, f (8) = O and ai = B+a, 
af > 


a, = 6% + a1 B + az. Conversely, suppose Re a, = 1m ag = Reag =Im aq = 0 
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and there exist x and 8 such that Re « = Re 8 = 0, 
2 


f («) = f (8) = 0 satisfying 322 + 2a, a -a + 4a, >0 
where 

a, = 8 + a1,a, = B2 + 218 + ag. 
Let 

K (A) = —, = A3 + a, A2 + a, A+ a, . SoK (a) = 0, 
and 

Rea, = Re B + Rea; = 0, Im a’, = Im 8? + a1 8 + az) = 0 


Rea, = Re (B2 + a; B2 + ap 8 + a3) = 0. Now case 3o0f Theorem 


3.2 implies that all roots of K (A) must have zero real parts. That ends the Proof. 


Finally, it is worth mentioning for the sake of applications that the arguments 
in all three theorems clarify the relations between the nature of the eigenvalues of the 
system and the coefficients of its characteristic equation. 
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A graphoidal cover of a graph G = (V, £) is a collection 4 of (not neces- 
sarily open) paths in G such that (a) every path in ¥ has at least two vertices 
(b) every vertex of Gis an internal vertex of at most one path in ¥, and (c) 
every edge of G is in some path in ¥. The graphoidal covering number y (G) 
of G is defined to be the minimum cardinality of a graphoidal cover of G. 
In this paper we determine the graphoidal covering numbers of trees, com- 
plete bipartite graphs, Hamiltcnian graphs and regular graphs. 


1. INTRODUCTION 


By a graph we mean a finite undirected graph without loops or multiple edges. 
We follow the notation and terminology of Harary2. All graphs considered in this 
paper are assumed to be connected graphs without isolated points. 


Let G = (V, E) bea graph. We denote the number of vertices in G by p and 
the number of edgesinG by g. If P = (uo, 1, uo, ..+, Un) is a path, not necessarily 
open, in G then uo and wn are called terminals of P and U1, Wo, ..., Un-1 are called 
internal vertices of P. We denote by ¢ (P) the number of internal vertices of P. The 
following definition of graphoidal covering number of G is given in Devadas Acharya}, 


Definition—Let G = (V, E) bea graph. A graphoidal cover of Gis a set ~ of 
(not necessarily open) paths in G satisfying the following conditions, 


(1) Every path in J has atleast two vertices. 
(2) Every vertex of G is an internal vertex of atmost one path in v. 
(3) Every edges of G is in some path in 

Let ( (G) denote the set of all graphoidal covers of G. Then y (G) = 


se | | is 
called the graphoidal covering number of G. Eg) 
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Thus y (G) is the minimum number of internally disjoint paths covering all the edges 
of G. We may further assume that every edge of G isin exactly one path of a 
graphoidal cover )sothatg = 2 | E(P)|=|v|+ 2 ¢(P). 

PEy PE 


In this paper we obtain the graphoidal covering numbers of trees, complete 
bipartite graphs, Hamiltonian graphs and regular graphs. 


MAIN RESULTS 
Theorem 1—y(G) = | E (G) | iff G = Ke. 


Proor: If G = Ke trivially y (G) = | E(G)| = 1. Suppose GA Ke. Let P 
be a path in G such that| E(P)| > 1. Then} = {P} U [E(G) — E(P)]} is a 
graphoidal cover of Gand || < | E(G) | so that y(G) < | E(G)|. 


Theorem 2—Let G be a tree with n vertices of degree 1. Then y (G) = n — |, 
Proor : We prove the result by induction on». Whenn = 2, G isa path and 


hence 7 (G) = 1. Suppose the result is true for any tree with n — | vertices of degree 
1. Let G be a tree with 1 vertices of degree 1 where n > |. 


Let P = (vo, v1, .--, ¥*) be a path in G such that d (vo) > 2, d(vk) = | and d(v) = 2 
fori =1,2,...,.4 —1. 


Then G] = G — (11, Ve, ..., ve} is a tree with n — 1 vertices of degree 1. Hence there 
exists a graphoidal cover | of G; such that | ¥| =n — 2. Clearly | U {P}isa 
graphoidal cover for G and hence y (7) <n — 1. 


Now suppose Pi, Pe, ..., Pn-2isa graphoidal cover of G. Thenp—1| = | E(G) | 
n-2 n.2 
= (n—2)+ = ¢(Pt). However 2 t(Pi)Qp—an which gives a contradiction. 
i=1 i=1 
Hence y (G) = n — |. This completes the induction and the proof. 


Corollary 1\—For any tree G, y (G) > A — 1 where A is the maximum degree 
of a vertex in G. 


Corollary 2—Let G be a tree with A> 2. Letv be a vertex in G such that 
d(v) =A. Then y(G) = 4 — | if and only if d(w) = 1 or 2 for all vertices w = v. 


Theorem 3—\Let t be a graphoidal cover of G such that every vertex v with 
d(v) > | is an internal vertex ofa pathiny. Then | ¥ | = y. 
Proor : Let be any graphoidal cover of G. Then Serie) ae. at 
PEy1 PEey 
Hence g —|¥1 | <7 —1 4 | so that | ti | 2 [d|. Hence| >| =». 
Corollary—\f there exists a graphoidal cover v of G such that every vertex of G 
is an internal vertex of a path in} then y (G) = q — P. 
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Pror: g= 2 |E(P)l=|¥l+ = t(P)=14l +P. 

PEy PEy 
Hence 
y(G)=|vl=q—-P. 

Theorem 4—Let G be a hamiltonian graph. Suppose there is a vertex v in G 
such that d(v) > 3. Then y (G) = q — p. 

Proor : Let C =(v = vo, 1, ---, ¥p-1, ¥) be a Hamilton cycle in G. Since 
d(v) > 3, there exist vertices, say x and y, distinct from v1 and vp-1 adjacent to ¥. 
Let P = (x, v, y). Let S denote the set of all edges of G not covered by C and P. 
Then v = {C, P} U Sis a graphoidal cover for G and every vertex of G is an internal 
vertex of apathin ¥. Hence y(G) = q — p. 


Corollary 1—For the wheel Wn = Ky + Cn-1 where n > 5, we have 
y (Wn) = q—p=n— 2. 
Corollary 2—For the complete bipartite graph Kn,n with n > 3, we have 
 (Kn,n) = q — p =n? — 2a. 
Corollary 3—Forn > 4, y (Kn) = q — p = 3 n(n — 3). 
Theorem 5—Let G be a graph with p > 5. If G has a Hamilton path P = (11, vo, 
..., Yp) Such that v; and vp have degrees > 3 in G, then y (G) = g — p. 
Proor : Similar to that of Theorem 4. 
We now determine y (G) for regular graphs. 
Any }-regular graph is Ke and 2-regular graph is a cycle. 
Hence y (G) = | if G is 1-regular or 2-regular. 
Theorem 6—Let G be a k-regular graph with k > 3. Then y (G) = q—p. 


Proor ; It is enough to construct a collection of mutually edge disjoint and 


internally disjoint paths in ( such that every vertex of G is an internal vertex of a path 
in the collection. 


Let Pi = (1, uo, ..., un) be a longest path in G so that all vertices adjacent to uz Or un 
are already in FP}, Since k > 3 we can find vertices x, y, z, w in P; such that x, y are 
distinct vertices each different from ue and are adjacent to v1 and z, w are distinct 
vertices each different from un-1 and are adjacent toun. If x, y, z, w, uy and un are 
all distinct, let Po = (x, u1, y) and Ps = (Z, un, w). If one of x, y coincides with un 


and one of z,w coincides with wu, say X = un andz = wy, let Me = (y, U1, Un, w). 


Thus we obtain a collection {P1, Po, P3} or {Py, P 


>} Of edge disjoint and internally 
disjoint paths in which u4, wo, .. 


» Un are internal vertices of one path in the collection. 
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If these vertices exhaust all the vertices in G, the proof is complete. If not let wi bea 
vertex not lying on P; and let Q; be a longest path in G containing w1 and internally 
disjoint with the paths constructed above. Ifthe end points of Qi are not in P; we 
make them internal vertices of some path as before. We continue this process until 
all the vertices are exhausted and we obtain a required collection of paths. 


We now proceed to determine the graphoidal covering number of complete bipar- 
tite graphs. Ki,1 is a path and Kg,9 is a cycle. Hence y (Ki,1) = y (K2,2) = 1. It follows 
from Theorem 2 that forn > 1, y (Ki,n) =n — 1. Also y (Ke,3) = 2. 


Theorem 7—Let G be the complete bipartite graph Km,n with m > 2 andn > 2 
orm =2andn> 3. Then y(G) = q -— Pp. 


PRooF : Case i—m = 2 andn > 3. 


Let ¥ ={v1, vo} and Y = {11, wo, ..., wn} be a bipartition of G. Let Py = (11,1,¥2,W2,"1), 
Po = (ve, W3, V1, Wa, Vo) and Pi-2 = (V1, We, V2) fori = 5, 6, ..., 7. 


Then | = {Pi, Po, ..., Pn-2} is a graphoidal cover of G and every vertex is an internal 
vertex of a path in ¥. 


Case ii—m = nandn > 2. 
The result follows from Theorem 4. 
Case iiicm = n+ 1 andn > 2. 
Let X = {v1, vo, ..., 'n,i1} and Y = {w4, Wa, ..., wn} be a bipartition of G. 
Let Py = (v1, W1, V2, We, ..., Yn, Wn, Vn,1), 
Po = (Wo, Vng1, W3) and P3 = (We, v1, W3). 
Every vertex of G is an internal vertex of one of these paths. 
Case iv—m >n+2andn > 2. 


Let X = {v1, vo, .... vm} and Y= {w1, Wa, ..., Wn} be a bipartition of G. Let 
Pi = (11, W1, V2, Wa, «++, Vn, Wn, V1), Po = (We, v1, W3) and Qt = (We, Vast, wg) where 
i= 1,2,...,.m —n. Every vertex is an internal vertex Of one of these paths. 


Hence y (G) = q — P. 
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The purpose of the present paper is to study some properties of «-Hausdorff 

subsets, almost closed mappings and almost upper semicontinuous decom- 

position. 

1. PRELIMINARIES 
Our notation is standard. No separation properties are assumed for spaces 
unless explicitly stated. 
A subset A of a space X is ‘regularly open’ iff IntCIA=A. A subset A of a space X 
is ‘regularly closed’ if Int CIA = A (Singal and Singal!'). 
A subset A of a space X is «-paracompact (“-nearly paracompact) iff for every open 
(regularly open) cover Q/ of A there exists an open X-locally finite family D which 
ee 
refines Q/ and covers A (Kovacevic® and Wine’?). 
A subset A of a space XY is «-paracompact (a-nearly paracompact) with respect to a 
subset & of X iff for every. open (regularly open) cover Q/={Ui : i EJ} of A there 
exists an open family WY = {Vj :j © J} such that: 
(i) D refines W, 
(ii) A CU{V;:j € J}, 
(iii) Y is locally finite at each point x € B. 
Subsets 4 and B of a space X are mutually “paracompact (mutually «-nearly 

paracompact) iff the subset A is «paracompact («-nearly paracompact) with respect to 


the subset B and the eye: Bis «paracompact («-nearly paracompact) with respect to 
the subset A (Kovacevic’*). 


A subset A of a space XY is «-Hausdorff iff any two points a and b, such thata € A 
and b € XA, can be strongly separated by open sets (Kovacevic). 


A subset A of space X is «-regular iff for any pointa € A and any open subset U of 


X containing a there exists an open subset V such that a € V coher (Kovacevic'). 
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A subset A of a space X is Lindeléf (nearly Lindeléf) iffevery open (regularly open) 
Vv 


cover ot A has a countable subcover Kelly’ and Kovacevic. 


A point p of a space Xisa P-point iff p € Int (1 (Un:n€ N}) whenever {Un : nEN} 
is a sequence of neighbourhoods of p (Kunnen’). 


A mapping f: ¥ > Y is almost closed (almost open) iff for any regularly closed (regu- 
larly open) set F of X, / (F) is closed (open) in Y (Singal and Singal''). 


A decomposition & of a space X is upper semicontinuous (almost upper semicontin- 

uous) iff for each D in D and each open (regularly open) set U containing D there 

exists an open set V such that D C V C U and Vis the union of members of & Kelly’ 
Vv 


and Kovacevic. 
2. ON «HAUSDORFF SUBSETS 


Definition 2.1—A subset A of a space X is «-almost paracompact in respect to 
a subset B of X iff for every open cover W%= (Ui: iE T} of A there exists an open 
family Y= {vi :j € J} such that: 


(i) Drefines QW, 
Gi) AC CI(U (V:7 € J)), 
(iii) Y is locally finite at each point x € B. 


Theorem 2.1—Let A be an «-Hausdorff «-almost paracompact subset with respect 
to each point of ¥\.A. Then A is closed. 


Proor: Let A be any a-Hausdorff «-almost paracompact subset with respect 
to each point of ¥\.A. Let a be any point of ¥\£4. For each x € A there exist 
open sets Ux and V, such that 


x € U,,a € Vz and Ux N Vz = 9. 
Now 
W = {Ux: x E A} 


is an open covering of A. Since A is u-almost paracompact with respect to each point 
of X \, A, then exists an open family Y = {Vi: j € J} such that: 


(i) Dis locally finite at a, 
ii) AC CI(U Wiss € J}), 
(iii) DW refines Y%. 


Since 7 is locally finite at a there exists an open neighbourhood U ofaand a finite 
subset Jo of J such that 


UNViF $b forj € Joand UM) Vi = ¢ forj © JNJo. 
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For each / € Jo there exists x (j) € A such that Vj C Ux(j). Let 
U, =U 11\(N {Vay : 7 € Jo}. 

U; be an open neighbourhood of a such that 
aE€uUucxXx\A 

hence, the subset Y \. A is open, i.e. the subset A is closed. 


Theorem 2.2—Let A be any «-Hausdorff nearly Lindeléf subset of aspace X. 
Leta € X \. A bea P-point. Then there are disjoint regularly open neighbourhoods 
of x and A, 


Consequently, if each point of X¥ \ Ais a P-point and Ais an a-Hausdorff 
nearly Lindelof subset of X, then A is closed. 


Proor: Since the subset A is «-Hausdorff, then for each point x € A there 
exist disjoint regularly open sets Uz and Vz such that 


x € Uz,a Ee Vz. 
Then 
YW = {Uz: x € A} 


is a regularly open covering of A. Then there existsa countable subset Ag of A such 
that 


AC U {Uz:x E Ap}. 
Let 
U= U (Uo: x € Az} and V = Int (-) (Ve: x € Ao}. 


Then U and V are open disjoint neighbourhoods of A and a respectively («(U) and 


a (V) — a(U) = Int CIU — are regularly open neighbourhoods of A and a respec- 
tively). 


Theorem 2.3—Let A and B be two disjoint «-Hausdorff nearly Lindeléf subsets 
of a space X such that each point of A and Bis a P-point. Then there exist disjoint 
regularly open neighbourhoods of A and B respectively. 


Proof: For each point x € B there exist disjoint regularly open subsets Uz 
and Vz such that 


x © Uzand AC Vz. 
The family 
QW = { Uz: x € B} 


is a regularly open covering of the subset B. Since B is nearly Lindelof then there 
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exists a countable subset By of B such that 
BC U Uz: x € Bo}. 
Let 
U=«(U {Uz: x € Bo}) and V = a(Int (N (V2: x E Bo }). 
Then U and V are disjoint regularly open neighbourhoods of B and A respectively. 


Theorem 2.4—Let A and Bbe any disjoint «-Hausdorff subsets of a space X 
such that : 


(i) each point of Bis a P-point, 

(ii) A is nearly Lindeldf, 

(iii) B is «-nearly paracompact with respect to the subset A. 

Then there are disjoint regularly open neighbourhoods of A and B respectively. 


Proor: For each point x € A there exist disjoint regularly open subsets Uz 
and Vz such that 


Vv ’ 
x € Uzand B C Vz (Theorem 2.3 in Kovacevic4). 
The family 
WW = {Uz:x € A} 


is a regularly open covering of A, hence there exists a countable subset Ao of A such 
that 


AC U{1Uz:x € Ao}. 
Let 
U = «(U {Uz: x © Ao}) and V = (Int (N (V2: x € Ao}). 


Then U and V are disjoint regularly open neighbourhoods of 4 and B respectively. 


3. Atmost CLOSED MAPPINGS 
v 4 . 
Using the similar method as in Kovacevic* we shall prove the following theorem. 


Theorem 3.1—Let X¥ be a topological space such that each point of Xis a 
P-point. Let f: X > Y be an almost closed mapping of the space X onto a space = 
such that : 


(i) for each point y € Y, /-} (y) is an x-Hausdorff subset of X, ii for each point yE Y, 
f-) (y) is nearly Lindel6f or «-nearly paracompact with respect to each subset f~! (z), 
z€ Yandz+y. 


Then Y is Hausdorff. 


Vv , 
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Proor: Let y1 and ye be any distinct points of Y. In any case, by hypothesis, 
there exist disjoint regularly open neighbourhoods U; and U2 of f~1 (yi) and f~1 (y2) 
respectively. Since the mapping fis almost closed, then there exist open subsets 
of Y V; and V2 containing y1 and ye respectively such that 

f-1 (yn) C f2 (1) C U1 5 f-1 (v2) C f-2 V2) C U2. 
Hence the result. 
Theorem 3.2—Let X be a topological space such that each point of X is a 


P-point. Let f: ¥ — Y bean almost closed mapping of a space X¥ onto a Lindeléf 
space Y such that: 


(i) for each point y € Y, f~1 (y) is an «-Hausdorff subset of X, 


(ii) for each point y € Y the subset f~! (y) is nearly Lindel6f or a-nearly para- 
compact with respect to each subset f-1(z),z € Yandz + y. 


Then / is continuous. 
ProoF: The proof is omitted. It is identical with the proof of Theorem 3.4 
Vv ’ 
in Kovacevic4. 


Theorem 3.3—If f is a closed and continuous mapping of a regular space X onto 


a space Y such that for each y € Y f-1 (y) is an «-paracompact subset with respect to 
the subset Y \, f-! (y), then Y is regular. 


Proor: Let y € Y and V be an open set containing y. Since the space X is 
regular and the subset f~1 ()) is «-paracompact with respect to the subset Y ts ed 6) 
Vv 


then, by Theorem 2.6 in Kovacevic4, there exists an open neighbourhood U of f~1 (y) 
such that 


EF (WiGiG SCI’ Caer, 
Since f is closed, there exists an open set W in Y such that 
y © Wand f-1(W) C U. 


Therefore, we have 


PGR ContO) tar Cl. ate 
Hence 


yEWCCWeY. 


Hence the result. 


Theorem 3.4—If fis an almost closed ma 
space Y such that for each point y € yf 


to the subset ¥ \, f-1 (y), then f is closed. 


pping of a regular space Y onto a 
(y) is an “-paracompact subset with respect 
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; Proor: Let A be any closed subset of XY and let y € X \. f(A). Since f~1 (y) 
C X \, A, there exists an open subset V such that 


FOO GEG CONG \.A. 
Then « (V) is a regularly open subset such that 
f(y) C«WV)CChaV)CX\A. 
Since f is almost closed, then there exists an open set W in Y such that 
ye Wandf3(y CfiW)Cca«V)CXN A. 
Therefore we have 
yEWCX\ f(A) 
hence Y \, f(A) is open. Then f(A) is closed. Hence fis closed. 


4. Acmost Upper SEMICONTINUOUS DECOMPOSITION 


Theorem 4.1—Let X be a topological space. Let D be an almost upper semi- 
continuous decomposition of X whose members are o-Hausdorff and mutually 
a-nearly paracompact subsets of ¥. Let D have a quotient topology. Then D isa 
Hausdorff space. 


Proor: The projection of the space X onto the space & is almost closed, 
a . 
hence, by Theorem 3.1 in Kovacevict. D is Hausdorff. 


Theorem 4.2—Let X be a topological space such that each point of Xis a P- 
point. Let D be an almost upper semicontinuous decomposition of X such that : 


(i) each member of & is an «-Hausdorff subset of X, 


(ii) each member DE Dis nearly Lindelof or «-nearly paracompact with res- 
pect to each-member V € QD, D=#YV, 


Let @ have the quotient topology. Then Q is Hausdorff. 


Proor: The projection of the space X onto the space D is almost closed, hence 
D is Hausdorff. 

Theorem 4.3— Let X be a regular topological space. Let D be an almost upper 
semicontinuous decomposition of the space X such that for each D € Q, Disan 
“-paracompact subset with respect to the subset X \.D. Let D have the quotient topo- 
logy. Then (i) D is upper semicontinuous, (ii) & is regular. 


Proor: The projection of the space X onto the space & is almost closed, hence 
the projection is closed. Hence D is an upper semicontinuous decomposition. Since 
the projection is closed then D is regular. 
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CONVEX UNIVALENT POLYNOMIALS 
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In this paper, we find a necessary condition and a sufficient condition tor a 
biquadratic polynomial to be convex univalent. 


1. INTRODUCTION 
Let f(z) = z + a2 2% + ... be analytic in the unit disc E = {z:|z| < 1}. Let 
F (z) be analytic and univalentin E. Let f (0) = F(0). If the image of the unit disc 
E under the mapping fis contained in the image of the disc E by the mapping F, then 
f (z) is said to be subordinate to F (z) and this is denoted by / (z) = F (2). 


f(z) = z+ a,2% + ... is said to be convex univalent in E if and only if 
Re {1 + 2f’ (z)[f (z)} > Oin E. 


f(z) =z+a,z2 +... is said to be starlike univalent in E if and only if 


Re {zf' (z)/f (z)} > O1n E. 


oo 
The Hadamard product or convolution of 2 power series f(z) =z+ 2 anz® 
n=2 


oo co 
and g(z)=7z+ 2 bn2"is defined as the power series (f * g)(z) =z # & anbn2". 
n=2 n=2 


Frankl’2 determined the necessary condition and sufficient condition for the 
polynomials oz + pd2 22 and oz +. pag 22 + Bag z3 to be convex univalent. The aim 
of this paper is to determine the necessary condition and sufficient condition for the 
biquadratic polynomial 5z + p42 72 + Baz z3 + yaq 24 to be convex univalent and to 
show that 


(12) 2 Va (2,f) 4 oz + pas 22 + oBas 28 + yaa 24 f(2) 
where V’3 (z, /) is the 3rd de la Vallee’ Poussin mean of the function f given by the 
formula V3 (z,f) = (3/10) z + (3/10) ag 22 + (1/20) ag z3. In general, the nth de la 
Vallee’ Poussin mean® of / is given by 


n 


Vn (2,f) =( 3 ) SS: ( Rots ) ax 2, (a, = 1). 


k=1 
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2. SUFFICIENT CONDITION 
Theorem 1—If f(z) = z + bz? + cz3 + dz4 where 6, c,d are real and non- 


negative, then the condition 


/ 


a2 


IV 


\ 


(2436); 


( 8c — 30d — 14cd ) 
4 


1+ Se — 32d 
ensures that f (z) is convex. 


Proor: It is enough to show that Re {1 + 2zf"(z)/f’(z)}} >0 in E. Put 
z=re?,O0OSr<1,059S32n. The above condition is equivalent to 


64d2 r6 + 84cd cos 6 r5 + (48bd cos 29 + 27 c?) r4 
+ (20d cos 36 + 30bc cos 6) r3 + (12¢ cos 29 + 852) r2 
+ 6bcos@r+1=0. eA ASD 


If r = 0, obviously (2.1) holds. Therefore assume that r #4 0. Denote the left hand 
side of (2.1) by F(@). Then 


F’ (0) = — rsin @ {84cd r4 + 192bd cos 6 r? + 30 (6d — 8dsin2 9 + bc) r2 
+ 48¢ cos @r + 6b}. 

F’ (0) = 0 gives sin 9 = 0 
or 

l4cd r4 — 32bd r3 + (30d + S5be) r2 — 8ecr +b = 0, we hana) 
We shall show that 

l4cd r4 — 32bd r3 + (30d + Sbe) r2 — 8er +b ei oa) 
Differentiate with respect to r and use the fact that d S 1/56. The worst case happens 


whenr = |. This is equivalent to l4cd — 32bd + 30d + Sbe — 8c +b=>0 and 
this gives 


» > Be — 30d — 14d 
= “T+ S5e— 32d * 


So when r lies in the open interval (0, 1), the inequality (2, 3) is satisfied. 


It remains to show that when F (6) takes its minimum value that (2.1) is satisfied. 
The value of F (6), when 9 = 0, is seen to contain only positive terms for r > 0, 


Hence we have to consider (2.1) only when @ takes the value x. When 6 = 7, (2.1) 
reduces to Pa 


(1 — 4br + 9er2 — 16dr3) (1 — 2br + 3¢r2 — 4dr3) > 0, »«(2,4) 
2, 
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The condition 


en 8c — 30d — 1l4cd 
it 1 + Se — 32d 


implies that b => (3/4) (8c — 30d + 14cd). This, in turn, shows that 1 — 4br + 9cr? 
— 16dr3 < | — 2br + 3cr2 — 4d3,0 <r <1. It is enough to show that 1 — 4br 
+ 9cr2 — 16dr3 has no root in the interval (0, 1]. Differentiate with respect to r and 
use b > (3/4) (8c — 30d + 14cd). This shows that the worst case happens when 
r = 1. This condition gives 1 => 4b — 9¢ + 16d. Therefore b < (1 + 9e — 16d)/4. 
Hence the theorem is proved. 


Corollary 1—Let 


3+ 14y am) 3+ 140 v one 1+ 90 y 


Py (2) = 4 — 10 20 


z3 aa vy 74 


and 0 S y $1/70. Then P, (z) is convex in E. 





F : Peas (2) *", 
Proor : P4 (z) is convex if Pete ree is convex. That is 
2(3 + 140 vy), (1+ 907) |, 4y P 
[wosiaiay ? SG+4y sey 
is convex. 
The equalities 
4y¥ | 
0s Sy =. 36 
and 
Ree fled =» 56.C1i 90 y) 
2(3 + 140) > 5: (3 +5 14'y) 34+ 14y¥ 5(3 + 14 y)? 
5(3+ 14y) — ; Foote the 200 128 y 





34+ 14y 3+ 147 
are equivalent to0 S y S (1/70). 


Corollary 2—Let f (z) =z + = an z™ be convex in E. Then so is 
n=2 


(Pa f)@) = Ott 


34+ 1407 5, 1+907 
in eR 2, 20 


X ag z3 + ¥_ a4 24. 


This follows using the celebrated result of Ruscheweyh et al.4. 
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eo . 
Corollary 3—Let f(z) = z+ 2 anz™ beconvexin £. Then 
n=2 


(Pn* f)(z) = Van(z,f) = f(z) 
and 
(1/2) z <x Vg (z,f) -< (Pa * f) (2) x Va (z,f) < fF (2). 

Proor : The relation (Pn * f)(z) ~ Vn(z, f) follows from Frank? and 
Ruscheweyh ef al.4 by taking f = Pn and 9 =fand b = Vn ( Zz. is ) To 
obtain the second reation, we note that (Py * f)(z) x Va(z,f) follows from the 
above relation forn = 4. Vs(z,f) < (Pa* f)(z) follows from the fact that 
(Pa * f) (z) is convex and thus in Theorem 6 of Frank’s paper*, we take 9 = Pa, 
~ = z/(1 — z)and f = V(z,2z/(1 — z). That (1/2)z-~ V3 (z,f) follows from the 
fact that the image of the unit disc by a convex function covers the disc of radius 1/2 
in the w-plane. 


3. NECESSARY CONDITIONS 


oo 
Theorem 2—\f for all convex f(z) = z+ Z an z™, oz + pag z2 + Bagz3 
n=2 


+ yaq z4 is convex and 

(1/2) z Vs (f, z) x92 + pag 22 + Bag z3 + vag z4 = f (z) 
then 

o = (1/2) + »—B+ yx and » S (1/6) + (8/3) 8 — (17/3) x. 


Proor : With f(z) = z/(1 — z) = z+ z2 + 23+ ... wecan find the minimum 
value of oz + pag z? + Bag z* + yay z4 on | z| = 1 whence | oz + pz2 + B28 + 
yz4;|>ao—p+68-y, If 


(1/2) z-< oz + pz? + Bz3 + yz4, we must have o — w+ B6—y S(1/2). With the 
same f(z) = z/(1 - z), if oz + pz2 + gz3 + yz4 - f(z) for real x, -l<xe< ¥ 
we must have ox + px? + Bx8 + yx4 > —1/2. Allowxto tend to — 1. We get 
¢S(1/2)+4—6+ y. Thuse = (1/2) + p— 8 + y. Forreala. -l<ca<l, 


consider 
1 + «2 — 
F (z) ={~* | 7 (a) - r (= 2 )] 
Sai Oe 1 — az 


where 7’ (z) = arctan z. Then F(z) is convex and maps F into the unit strip 


1+ 22 


1 — «2 


+ a2 
ial 





{arc tan « — w/4} S Rez < 





{arc tan « + 2/4}, 


The function F(z) has the Taylor expansion 
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oo 





F(z)=z+ SS BO ee aye 
n sin @ 
n=2 
where 
c 1 — «2 2a 
on ak and cos @ = ice 


If the polynomial {(1/2) + » — 8 + y}z + paz z2 + Bag 23 + yay z4 ~ F(z), we 
must have for real x, —1 < x < 1 


arctan a — 7/4 S {1/2 + » — B+ y}sin 0x + (p/2) sin 26 x? 
+ (8/3) sin 39 x® + (y/4) sin 46 x. 
But arc tan « — 7/4 = — 6/2. Taking the limit as x tends to —1, 
—@/2 S —{1/2 + » — B + y} sin @ + (4/2) sin 26 — (8/3) sin 
36 + (y/4) sin 40. 


Taking the limit as « tends to 1 corresponds to @ tending to zero. By repeated applica- 
tion of L’Hospital’s rule, result follows. 


Theorem 3—TIf for all convex f(z) = z + > an z” 
n=2 


3 + ly 
4 


3 + 140y 


as 10 


a2 z2 + a 


0 a3 z3> + y aq 24 


is convex and 


esti SFA gas ot ASO 


4 10 
+ yagz4~ f(z) 


(1/2) z = V3(f, z) R 


then y S 1/70. 
Proor : Since os = 1/2 + »— 8+ y when p = 1/6 + (8/3) 8 — (17/3) y, the 
polynomial f(z) = oz + paz z? + Bag z3 + yag 24 becomes 


eB tage ee ol at 1 + 90y 


2 
4 10 a2 2° + ——79 


a3 z3 + yaq 24. 


The required result is a special case of Theorem 6 in Frank’s paper? for n= 4. The 
proof is complete. 
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(In Memory of L. Shri L. S. Bosanquet) 


The object of the present paper is to determine a class of sequnces 4 = (An) 
and absolute summability methods { A | for which | A| and its modified 
methods | A’, | are equivalent. The special case of Noérlund matrices is 
taken for special study. 


1. INTRODUCTION AND NOTATIONS 
Let c and 1 denote the set of convergent and absolutely convergent sequences 
x= (xn). Given an infinite series z an With (sn) as the sequence of nth partial sums 
and an infinite matrix A = (ank), we write the sequence-to-sequence transformation 
= A(s) as 


Cc 
tn = XZ nk Sk oe ISH 
k=0 


assuming that tn exists foreachn > 0. We write On = tn — tn-1, t-1 = 0. The 


series 55 an is said to be summable (A) to the value sif ¢ € cand z On = s; and 
n=0 n= 


is said to be absolutely summable 4 if ® € 1 i.e. x | On | < co, (Summation with- 
out limits means summation from 0 to 9). 


For a given sequence A = (An) let tn be the sequence of A-transformation of 
(An an), that is, 


™ = = nk Ak ak. et f Bey 


Suppose that An + 0 whenever n > 0 and let 


B = (bnk) =( “ae ) (1.3) 





in the case Ay = 0, we define ao9 = boo. We write 


oo 
vn = tn/An = LZ bnk ak. w--(1.4) 
k=0 
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The series an is said to be summable by the ‘modified A-means of weight’ A or sum- 
mable (4’, A) to the value sif 2 dn = s; and absolutely summable by modified A 


means or summable | 4’,A | ify € 1. In the case 





(1.5) 


...(1.6) 


os( ded) 


.. (1.8) 


x= On = Tn 
we say that the method A and | A | are identical to their modified means with weight 
A. 

If the matrix is lower triangular, i. e. ank = 0(k > n), then (1.1) can be putin 
the form 
n — 
tn = & ank ak 
k=0 

where 

= n 

aAnk = = any 

Pak 
and 
nea 
On = tn — tn-1 = Z ank ak 
k=0 

where 

A A A = = 

A = (ank) and ank = @nk — an-1, k. 
Thus in the special case when 

. Ak ank 

ank = —— =, Dak 

An 


the method (A) is identical with (A’, A) mean. 


Let (A; ) be the sequence of Cesaro coefficients determined by 


Z Ae x® = (1 — x1 | x] < 1). 
If 





(kK <n) 


0 (k > n) 


and « > — 1, An =n, then (1.9) holds; for we know that 


(1.9) 


a, OS UN) 
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A~ ‘ 
n-k “Denar k 
ee er gt Ce gee 
a a, nn An-k : this 
n n-1 


Thus the Cesaro method (C, «) is identical with its modified mean. 
The method (A) is said to be ‘equivalent’ to its modified mean with weight A if 
x Dn converges <> X Un converges. os oO ea 
Similarly | A | is said to be equivalent to | 4’, A | if 
PElSeVE l. ; ...(1.13) 


When (1.12) [respectively (1.13)] holds we write (A) ~ (4’, A), 
[respectively | A | ~ | 4’,A | ]. We write 


@ (A) = {AE R: (A) ~ (A, Aj}. ...(1.14) 
Let N be a Nérlund martix (JN, p) defined by 


a (k <n) 





Pn 
0 (k>n) PPCHTS) 


Onk = 


where Pn = po + pi +... + Pn #O forn > 0. Let (N’, p, A) denote the modified 
Norlund mean. 


When p = (pn) satisfies Kaluza condition : 


pn > 0, Pnyi pani > P, , Pnyi < Pn (1.16) 


we write p © GH. Das‘ has shown that when p € DX, 
nE€ 6(|N{)i.e. |N, p| ~ 1’, p,A | with An = 7. Pe LEY 


Bosanquet and Das! have shown that a Norlund method is identical wth its modified 
Nérlund mean of weight An = n if and only if the Ndérlund matrix is a 
Cesaro matrix. 


It is easy to verify that Riesz method (R, pn, 1) is identical to its modified mean 


with weight 


nN 
An = r 
Hn41 — Bn 





When a method is identical with its modified mean as in the case of Cesaro or 
Hausdorff mean [see Hardy®, p. 247] it is an ideal situation; because the modified 
means are usually simpler to work with. Even when identification fails, equivalence 
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may sometimes hold and this serves as well. There are again Noérlund methods which 
are not equivalent to their modified means!. 


The knowledge of the scope of @ (A) may be helpful in dealing with the summa- 
bility factor problem: Suppose that we are required to determine if © en an is summable 
A(or|A|). If we know thataA € @ (A), then instead of considering A-mean of & en 


an, we may consider 
] 
(an) = ( — 2% dank €k Ak ar), 
An 


Now we may choose (en), or (An) € 0 (A), Or both, in such a way that 
en An = nor én An = fF 


In that case «n reduces respectively to 


én 
a’ = —  ank kar 
n n 


and 


a” == €n & Gnk ak 
and obviously ~’ , «* are simpler to work with as the factor én is outside the sigmas. 


We need some additional definitions and notations. 


The matrix A = (ank) is said to be normal if it is lower triangular with non- 


zero diagonal elements. It may be noted that the normal matrix A hasa two sided 


normal inverse which is denoted as A~! = (a; 


We write 
™m m 
a -1 * -1 
m= > at, of, = ealgcs wea(1.18) 
k=v k=v 
co oo 
= -1 ers =] 
Oy ss CP 6 = SS | a, | Be Pal) 
k=v k=v s 
whenever these exist. 
We write 
Pn-k Pn_k-1 
» fork <n 
Q (n, k) oe Pn Pa ee 


0 ,fork>n. -.(1,20) 
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Let the sequence (cn) be defined formally by 


(Z cn x”) (Z pn x”) = 1 set i2l) 
where po ~ 0, that is, by equations 
n 
oy ays ee 1 (a = 0) (1.22) 
oe 0 (n > 0) 
and let | 
c_4= 0. 


We write, for any sequence ( fn) 


A° fn = fn, A fn = fn — fay, A” fn =O (A*-1 fn) 

Vo fn = fn, V fu = fn — frat, Vr fn = VV fr) 
forh = 1, 2,3,..., and define y! fn =fo +fit... + fn. We adopt the con- 
vention that fn = Oifm <0. Further we write 


h h_ h- as 
f= = f= Sera ae as -1) me nie 1) 


for h = 1, 2, 


In Section 2 we make an attempt to discuss about the general mean | 4’,A|, 
whereas in Section 3 and onward we confine ourselves only to modified Norlund mean 
| N’,p,A|. The question of modified (A’, A) mean has however not been discussed 


in the paper. 
2. Matrix METHODS 
We first obtain the following fundamental lemma. 
Lemma 1—Let A = (ank) be a triangular matrix and let A and B be defined by 
(1.3) and (1.8). Suppose also that A 1 and B-1, the inverses of ‘ and B respectively, 


exist. Let 


n 


4 B-l = E> dnk bt ) = G = (gn) 


k=r 


and 


n 
B A} =( >, bnk PB ) = H = (hn,,). 
ker 
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(a) In order that | A’,A | > | A | itis necessary and sufficient that, for all ie 


S lg laek (2.1) 


n=k 
where K is an absolute positive constant not necessarily the same at each occurrence. 


(b) In order that | A | > | A’, A |, it is necessary and sufficient that, for all k, 


8 


| hn | < K. 


ll 
> 


n 


Proor : We write (1.4) and (1.7) as 
v= Bla), © =A (a). (2.2) 
Since the inverses B-!, A-! exist, we obtain 


a = B-1(¥), a = A-1 (0), 


Hence 


® =-AB1() = GW) .a-(2.3) 


and 


Y = BA1() = H(9). ...(2.4) 
It follows from a result of Knopp and Lorentz® and (2.3) that 
vEl=OE1(i.e.G:1—> 1) 


if and only if (2.1) holds. Similarly (2.4) yields ® € 1 = ~ € lif and only if (2.2) 
holds. This completes the proof. 


Now we prove 


Theorem 1—Let A bea normal matrix and (A) be a positive and non-decreasing 
sequence if 





AAT? [os +23) 

Ae (2.6) 
An — Xk +. Akii — Ak 

DS 7 | Ak ank | = o( FST ONG Mi <2 Gea 

n>k 


and 
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S lone |4 ( 5-)=0( 5 ) ..(2.8) 


where (@nk) is defined in (1.18). Then 
[a AD pot fb. 


Remarks : (i) Note that the condition (2.5) means that if x € 1, then A Ama x 
€ 1, and then this is so if and only if 


co n 
SS | > any a) | < K, for » = 0, 1,2, ... (See Knopp and Lorentz‘). 
n= v= 


.. (2.9) 
Similarly (2.6) holds if and only if 
an A 
DY | anu | < K wank det} 
n=" 


for wo 0, 1,72, 32. 


(ii) In Lemma 3, we have discussed some simple conditions under which (2.8) 
holds good. 


Proor : To prove the theorem, we have to prove (by Lemma | (a)) that, for 
bet OS 15 2e 223 


n 


ioe) co , 
Jn => Lim t= >| > anv bgt 





[Penk 
n= n= {=F 
Since 
255 | Au -1 
6 oom Av ap 
Au -1 (An = Av) Au ] 
See Rage SSS 
we obtain 


TP ime bagi a bo 


where, by hypothesis (2.5) and since An > 0 and non-decreasing 


[ova) Co 
1) cag mee or 
x eri A > hg | any ay | 
ann ver 


(equation continued on p. 354) 
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a A 
<a S any a 1<K 
n= 





00 n 
1 at. Pee 
IP m= de Dap | Dam ae | 
n=- v= 
co n 
1 Ay — Av * 
= Die | Do (AG ae ) Oe I 
n= v=P 
oo \ n ‘ 
< mu > oe, | SS An A (1/Ay) anv41 Avy | 
n=p y= 


co n 
l An—A A 
Tee > ao >*— * (Av dav) Bop | 
n=-~ v= 
= the + Aad say. 


Now by (2.6) and (2.8) we obtain 


co 1 co 
‘hee OS | @vu | A ( x) SS | anv 
v= , n=v 


co 


<K Ae > bow 1 a(s ) 
v 


v= 





= K,for p = 0,1, 2... 
and by (2.7) and (2.8) 


ioe) 
uae Mm > | Avy | 


ia 
Mi 
rs 
~ 
S 
a 
=> 





i 
y= n= 
co 
<K u> dane a ( : )<kx 
< se 


for» = 0,1, 2, .... 
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This completes the proof of Theorem 1. 


Remark : It follows from Lemma | (b) that a necessary and sufficient condition 
for | A | > | A’, A | is that (2.2) should hold, ive. 


Clearly (2.2) holds if 
co 
Gy) 3 | bas | = K,'k = 0, I, 2, ... 
n=k 
and 


oo 


A 
Gi) > an 1k =G 1,2... 


v= 


However the above condition (ii) does not hold in the Cesaro case. For in that case 


n 
A 
=] % —%—2 
a —_ : A, = aed : 
r=k 


It may be remarked that so far it has not been possible for the author to determine 
some suitable simple conditions on a general matrix A in order that the condition (2.2) 
for | A | + | 4’,A | be satisfied so as to include both the Cesaro and the general 
Nérlund case. However the problem of establishing | N’, p, A | ~ |, p| has been 
tackled in a somewhat satisfactory manner in later sections. 


An analysis of Conditions 


We now make an analysis of conditions in Theorem |. For this we make some 
preparation. 


Let 


We write A € D7* if 


Ais normal, ank > 0 fork <n AeA) 
dnk < dnyke-1for0 <k <a, bon 200z) 
dnk <1 forO <k <n. bs (oel3) 


The following results are known. 
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Lemma 2—(a) Let A satisfy conditions (2.11) and (2.12). Then its normal in- 
verse satisfies 


=] aL l 
nk <O0(k <n), Cage ie 


> 0. Wale} 





(b) If A satisfies conditions (2.11) and (2 13), then 
n 


>, ob 20,0<ren,n=0,2,... (2.15) 


ker 


(c) Let A € DP *. Then 


n 
> | a | < 2Jann .. (2.16) 
k=0 


n 
(d) Ifin addition, = ank is non-decreasing as n increases, then 
k=0 


fw ofa age cel Ge ee eT ...(2.17) 


See Peyerimhoff® (p. 33) for the result (2.14). For the result (2.15), see Peyerimhoff? 
(Satz4), The result (2.18) isa trivial deduction from (2.14) and (2.15). The result 
(2.17) is due to Russell9, 


Remark : The result (2.15) fails to hold under the hypotheses (2.11) and (2.12). 
For example, let (an) be a sequence of positive numbers and define 


aAnk = Qn (0 <k < n). 


Then the hypotheses (2.11) and (2.12) are satisfied. But if (tm) is the A-transform of 
(sn), then 


tn = Qn ZX ak 
k=0 
so that 
sa=-> In-1 tn x 
an-1 an 
Thus 
if l 
Fae (k =n —1) 
bor } | 
nk | on (k cee n) 


| 
| 0 (otherwise). 
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Hence for a given n > 1, (2.15) holds for all k < nif and only if 
an-1 = an, 
In the following theorem we examine the case for dnx a constant. 
Theorem 2—Let A be normal and let ank # O (k < n) 
Then 
dnk = «(4 0), (k <n) ... (2.18) 
if and only if 
a) =O(n>v+ 2) ! 
\ 2.19) 
Ak 41, k41 ee — : 
On the otherhand 
a) =O(k+25n<N+ ) | 
and . ...(2.20) 
-1 
Goeae 7? j 
if and only if 
dnk =a(k <n <N) 
date nae ect ayi: (2.21) 
Proor : Suppose that (2.19) holds. Since 
at, REF Bia yy rik = O ...(2.22) 
we have 
dik = Sede ee 


-1 
— Qk4,1,k41 @ 4 = &. 
akk +49 k+1’k 


Also we have (see Russell®, p. 101) 


n 
= ait 
An41,%41 aoe = > any (dnk — dnv) a). - ...(2.23) 
v=k+1 


Since eos = O(n>k+}), (2.23) reduces to : 


pits dar) Saae O FN) (2.24) 


and since forn > k +1, 4n,k+1 = Beccie + 0, we obtain from (2.24) that 
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dnk = dn,k41 (n Sm k +1). 


Hence (2.18) holds. 
If, on the contrary, (2.18) holds, then 
Ank = atk agr 
and from the formula (2.22) and (2.23) we can reduce (2.19). 


If (2.20) holds, we obtain 
ON +2,N42 ayes = 4n41,N41 (dw41,4 — dn41,k41) Pe 
and from this we obtain (2.21). The proof that (2.21) implies (2.20) follows from (2.23). 


This completes the proof. 
We write A € D7 ** if A satisfies conditions (2.11), (2.12) and the following : 
Gnk S ankyi (OS Kk <n — 1), PAS eae RY 
Now we prove : 
Theorem 3—Let A © G**. Then for all m, n and ih ahaa § 


(i) Omn > Omyi,n > 0; On > 0; 


(ii) ann (6* + On) = 2; 


r 
(ili) 0 < SS ank a, S any Ory (0 Sv < rs n); 
k=v 


m 00 
; > Ss -1 
(iv) akn | a Ber 
k=n v=m+] 


Proor : Letr En. Since ank S Qn,ky1 and og SO0(kK>yv +1) by Lemma 
2 (a), we obtain: 


r 





as | 1 5 
> ank a) =a) anv — Were | kid ee a | nr 
k=v 
= a ee Anyv4 -1 a 
ow ( 0 — lal eappeaeeee tanr a 
S Any aw ni | a | -1 
vy veiw 1 — 14 |) = any ry. ...(2.26) 
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On the otherhand 


r n 
—1 ps -1 ] (n = v) 
ank a, => Qnk a, = Sny = 
> ky >: 7 Be ny ‘o (aes) 


Hence from (2.26) we obtain (iii). 

Now from (iii) we obtain @mn > 0 and by making m — oo, we find that @n 
exists and that @n > 0. Also @mn > 0m,1,» is trivially true after Lemma 2(a). 

Since for n > y, 


n n 
> -1 5 =] -1 
a, a Ory = 2a, a8 0 
k=v k=Vv 


we obtain, by making n > oo, 
2 


Hence (ii) follows. Again 





m ©o oo 
== oS ky ae my ke 
k=v k?=mvr1 k=m+1 
so that 
(ore) 
la | <e 
kv — my 
k=mt1 
Hence 
™ co ™ 
> akn > | at is > akn @mk = |. 
k=n r=mt1 k=n 


This completes the proof of Theorem 3. 


As commented earlier, we now give the following result before we take up the 
study of the problem for Nérlund means. 


Lemma 3—Let A € GPU ** and let An > 9, non-decreasing such that 


n (Any1 — An) = O (An). ..0(2.27) 
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Also let 


Spee OLY: (2.28) 


ver 
Then (2.8) holds. 
Proor : By Theorem 3 (i) 


Qu 


l 
O< O2p- 5p < a > Avy 


v= 


and therefore (2.28) implies that 


O2p,n = O (I). 
Now by Theorem 3 (i) 
co 2¢ co | 
l ete 
> (a )-(S+ > ) ma (a) 
k= kab k=2e+1 
2u co 
= Okp i 
= o( Tis ) + Gm» > a( 5 ) 
k=P K=2Qe+] 
i< ] 
] 
=0(5 ) r 2S Okey + Oay,u. wn (4. ) 


by hypotheses (2.27) and (2.28). 


This completes the proof. 


3. Mopirie>D NorLuNpD MBANS 
We shall now study conditions on A and p such that Np, Allee Noe 


Recall that 2 an is said to be summable pay <p. 
NSP, AY at 


Slee 
An 


A | and written as Dan € | 


where 
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We may define (N’, p, A) method to be absolutely conservative if 
Z| an| < 00 > Zan €|N',p,al. 


It is easily verified that the modified Nérlund means (N , p, A) is absolutely conservative 
if and only if 


= . | pak | 
Jk = Ak >-4 eet <a K tst) 


n=k 
for ki 0, 1,2, 2.. 
The following Lemma gives simpler conditions in which (3.1) holds. 
Lemma 4—Let (An) be positive and non-decreasing and such that 
00 
> a = o( x) G2) 
n=k 
and let (pn) satisfy the conditions : 


() PS = > | pe| = 0(| Pal) 


k=0 
(ii) (2 + 1) pn = O( | Pn] ). 


Then (N’, p, A) is absolutely conservative. 


Proor : From (3.1) 


\ Ak | pn-x | 
n-(S+ s* ea 





n=2Kk+1 
~ | pnt | 
K Pn-k 
< ay "eae Bae ae S An | Pn-z | 
n=k n=2k+1 
< 1 
. — — <K 
<k+Ku > (in — k) Anok 
n=2k+1 


by hypotheses. 
This completes the proof. 


Remark: The hypothesis (3.2) is automatically satisfied if (An) is positive and 


362 G. DAS 
there is a constant « > 0 such that for n > 1 


An \. 
Pe non-decreasing. 


For 





— 1 < a 1 
n®* n 
» Oy ee >, n™"1 An < >, n=t1 -o( Ak : 


n=k n=k n= 
We now obtain the following fundamental lemma. 
Lemma 5—(a) In order that 
|N’,p,A| > 1N, P| 


it is necessary and sufficient that, for k > 0, 


n 


0 
Jk =)Xx | Px | >! Sy O(n») cv-4 | SK 


nk u=k 
where Q (n, v) and (cn) are defined by (1.20) and (1.21) respectively. 
(b) In order that 
[N,p|>|N’, pA 


it is necessary and sufficient that 


co n kj 
ee S ae oh »S AsiRe? eS hela Lin 5 
n=k v=k 


h=y 
Gy Ree) 9 36 Se 
Proor: When A is the (N, p) matrix, we have 


tec { Pk cn-k (k < n) 
nk 0 (k > n); 


A 


ank = © (n, k); 


Ak Pn-k 
bak = — —? 
oe An Pn ” 
n 
Ay 
SS  , DS Py (en-+ — Cn-y-1) 
r=k 


K-1 
es x Py (cnr — Cn-r-1) 
r=0 


(equation continued on p. 363) 


MODIFIED MEANS 363 


k_1 
= Pe_ycn-k— & Ppren_s; 
r=0 


—< Ak Pk cn_k 
| Ae Galen Saal tetee lag 
nk an 


The lemma now follows from Lemma |. 


It may be noted that the condition pn € D7 covers the Cesaro case pn = ae 


0 < « <1 whereas the condition V” pn € DV covers the Cesaro cases: h <a <h 
«oe ee Sat Be ate 


Now we state our main theorems. 


Theorem 4—Let An > 0 and increasing to co. Leth bea non-negative integer 
and let 


yi pn © DT ...(3.3) 
and 


not rw net ( 5 ) = 0.0). (3.4) 
An 
In the case h = 0, we further assume that (3.2) holds, i. e. 


> = ot terete 


n=k 





Then| N’,p,A|>|Npl. 


Theorem 5--Let (An) satisfy (3.2), hold. Let A be a non-negative integer and 
let y* pn € DZ. Also let 


nh+1 Atl An = O (An). sat ais) 
In the case h = 0, we further assume that either 
> si _ ( cide’ ) .-.(3.6) 
nPn Pk 
n=k 
or 
( Ae) is non-decreasing. oe(3./) 


Then 
|\N,p|>IN,pAl- 
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Remarks: (i) It may be seen that a positive sequence A condition (3.2) is 
automatically satisfied whenever (3.7) holds. 


(ii) The hypotheis (3.6) is satisfied in the case Pn = A, ,« > O and is not 
n 


: l 
satisfied when Pn = x ser a a 
v=0 





Theorem 6—Let (An) satisfy (3.2) and let An > c0, n > cco, Let h bea non- 
negative integer such that ¥" pn € DV and (3.4). hold. Inthe caseh = 0, we fur- 
ther assume that either (3.6) or (3.7) holds. Then | N, p| = |N',p,A|. 


4. Lemmas 
We need the following lemmas for the proof of Theorem 4. 


Lemma 6—Let h be a positive integer, y* pn = gn and Qn = go+ git ... 
+ qn. If gn > 0 and 
Onjqn is non-decreasing (4.4) 
then 


Pnjpn is non-decreasing. ..-(4.2) 


Proof: It is enough to prove the lemma for h = }, the general result will then 


follow from the repeated applications of this case. Changing the nctation in an obvious 
way and writing 


Pic Ppt Party te Be 


it is enough to prove that if pn > 0 and Pn/pn is non-decreasing, then PP, is non- 
decreasing. The result to be proved can be put in the form that 


PO), Pa — Pax P® > 0, (4.3) 


We prove (4.3) by induction on n. It is trivial that (4.3) holds when nm = 0, Let 
now m > | and suppose that (4.3) holds with n replaced by nm — 1, 


We have by hypothesis 
Pn oa Pn41 
POR Rms anaes 
Pn Pny41 


so that 


Pn (Pni1 — Pn) < Pnyi (Pn — Pn_y); 
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that is to say 


Pay Pa <P. . (4.4) 


Thus by (4.4) 
PO Pa —Pnyr PO = (PO + Pays) Pa — Pas Py 


= P| Pp + Pair (Pn — P%) 


= PO Pa — Pay P™ 


n-1 
Pp? 
~<a i) re (1) 
a a Pn : er. 1 


by the induction hypothesis. 


This completes the proof of the lemma. 


Lemma 7—Let h be a non-negative integer and let A” pn € St. Then 


(i) oe == tq > 0, a <= 0(n > 1); 


(iy eO*) > cht SO forn > 0; 


co 
(iii) > a x" js absolutely convergent for | x ae Fs 


n=0 
lo a) 
h h 
(iv) > feoii<nce 3 
k=n*1 


(v) (Vv?! pn) ee <1; 


(vi) (V%1 pn) ae < 2n +1; 


r 


h h h (41) 
Pi OS > Gs ees SO, Pere OSES <n". 


b=k 


Proor : This follows from Kaluza’s theorem (see Hardy®, Theorem 22), the 
identity (obtained from (1.21)): 
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co 


(> Caan ) BS (v" pn) x") = 1 (4.5) 


n=0 n=0 
and from the Lemmas 3, 4, 6 of Das4 (see also Das*, Lemmas | and 2). 
Lemma 8—Let y* pn © DZ. LetO &s < handO <r<ch. Then 
(i) ys pn > 0, non-decreasing and 
V*> pn > 9° as Nn > O°; 
(i) (n+h—rt+1) 9 pn <(h —4 + 1) YT! po; 


ie . Al(n +h —r)! 
(iii) "pn <G+Diho nL pn 





(A+ Il(n+h—r)! Ps <0 ( Pn } 


~ (h-rln+th+ Dd! i ut 


n 

; h— = 

(iv) Pn = 5) Ay, (A pk) < (yl pn) A. ; 
k=0 


and 


n 
h h 
am >. A,-~ V' Pk > V" pn. A,” ; 





k=0 
(v) rie - 0( 2); 
vil pn. Fa j* 


(vi) The method (N, p) is regular and absolutely regular. 


ProoF: Since 
n 


_ h_-r1 
V" pn -> Nl, ..-(4.6) 
v=0 


P h-r- 
and since y” pn > 0, A, are 0, the result (i) follows. To prove (ii), we proceed as 


in Das* (Lemma 3). It follows from (4.6) that 


(h—r+ iIyt pn -(nt+h-r+1)¥" pn 


= > (h—-r+ DA. _@th—r4)) An) }y" py. ...(4.7) 


v=0 
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But 
n 
SS Gar) e fata repay 
v=0 
h_r+ 
—(hort at? -(@th—orti ae =0 (48) 
as 
h_-r+1 
An on thr tl 
rege h—-rt+i 


which is an increasing function of n. Hence it follows that there exists v <vo (”) 
such that the expression bracketed in (4.7) is non-negative for v < vo and negative for 
» > vo. Hence the right of (4.7) is greater than or equal to 


n 
xc —r+ Ny pee — (n+ h— F+ 1) Via } VA Pro 
v=0 
which is 0 by another application of (4.8). Now (ii) follows. (iii) follows by repeated 
application of (ii). (iv) is obvious. (vi) follows from (iii) and (iv). 


Since pn > 0, the necessary and sufficient condition for (N, p) to be regular is 
that 


pn = 0 (Pn) ...(4.9) 


and this holds because of (iii). To prove the absolute regularity, in view of (4.9) it is 
enough to show that, uniformly in vy 2 9, 


S| O(n») | = 0 (1). (4.10) 


n=v 


Since 
Pn-+ i Pn _v-1 
Pn Pn-1 





Q (n, v) = 


_ Pn Pn-vy — Pn_v Pn 


is Pn Pn_i 


( Pu - Pn_y ) 
Een pn Pav 
Pn Pn-1 


it follows that, as pn 7 


0. 
Q(n, v) > 0 
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if and only if Pn/pn is non-decreasing. But since y’ pn = gn is non-increasing, it 
follows that Qn/qn is non-decreasing. Hence, by Lemma 6, it follows that Pn/pn is 


non-decreasing. Hence © (n, v) > 0. So 
co 


co 
: Pm-_» 
Q (n, -S$ 20, = lim = | 
n= 


n=Vv 





by use of (4.9). 
This completes the proof. 


Lemma 9—(a) Let h be a non-negative integer such that y” pn € DT. Ifh = 0 
we further assume that (3.6) holds. Then 


oth 
Doe -pet * O (1). 
v=Kk 
(b) Let y* pn © DT and let (An) satisfy (3.2). Then 


Ak Pk ee Iyhtiq, > O (1). 


(241) 


Proor: Leth > 1. Since C., is non-increasing we obtain 


men 


u=k v=k v=2k+1 
(h,2) 
CK 


co 
; (241) ] 
Sk + net + % > ) (+ 1rd 
v=2k+1 





cs o( at I O l 1 
vin Jaro +°0(—ap )eaap 


= O(1]Px) 
by Lemma 7(v), (vi) and Lemma 8 (iv). 


if h = 0, then we split the sigma as before and obtain 


= hie ee 

—F <>— + Saks canlt 
Wt ~ k+1 = G+ 1) Pye 
v v=2k+1 


me C1) Fe) 
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by Lemma 7 (v), (vi) and hypothesis (3.6). 
This proves (a). For the proof of (b), we proceed as before and obtain 


oth) Aa 
2, fe ae 
v=2k+1 
l 
= o( De Pk | 
by (3.2). 


Lemma 10—Let A be a non-negative integer such that Vy" pn € DG. Let An > 9 
non-decreasing. Then 





06 Ai. Pn_k 


| l 
(i} > Ol aeqap ) for = ts 


n=k 


Ar Pn_-k 1 
(ii) Saas (n+1)Pn = 0( cr) for r= 0, [; fs oe h; 


co Ar Pn_-k 
An Pn 


n=k 


= 0( apr) fore = ie te oy 


(iii) (k++ 1)? Ak 


The result (iii) remains valid for r = 0 with the additional restriction contained in 


(3.2) [see Lemma 4]. 
Proor : By Lemma 8 (i), (iii) 


(ore) Ay. pn -k oo fg , Pn- k 
oe -(>+ 3) 
ak n=k n=2kt+1 
Pn-k k 
< a See pn-k + K (n—ky*l Pn 
n=k n=2k+1 
_ l 
vy" 1 pk ee es 
Nae Segue a (n—ky 
n=2k+1 
K 
< (A + 1) © 


The results (ii) and (iii) are immediate corol- 


; ts the proof of (i). 
his compres It (ii) can be proved as in Lemma 4 


laries of (i) forr = |. Whenr = 0, the resu 
by use of Lemma 8 (iii). 
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Lemma 11—Let / be a non-negative integer and let y” pn € DZ. Then 


Ba a, a(n,k) |= 0( a) 


+ 1)r 
fory. == 001-24 


Proor: Thecase r = 


0 is (4.10) which has been proved earlier in the proof 
of Lemma 8 (vi). 


Consider the case r > 1. 





Since 
: Me Pn-_k Ma Pn_k-1 
Aus ew a eee 
Bei ) Pn Pn1 
A’ pa_t 
= aaa ge —_ ___ps At} k 
Pn Pn Pn_1 beg Pare 
we obtain : 
ae ka Ai Pn-k oo 
SS JArTQ(n,k)|< > —_— + > = Ps wis a nn Bas Pn-k_1 
n=k n=k : 


n=k 


O ] — A pn-k_] 
= ( (k+l) as o> (n + 1) Pn 
n=k 
Pa o( Boe sents 
ct) 
by Lemma 8 (iii) and Lemma 10 (ii). 


Lemma 12—Let h bea non-negative integer, v4 pn € GP and let An>O be 
such that An increases to co and 


n (Anyi — An) = O (An). 








Then 
fore) hai 
gk= Dy (Anghsy — Ak,n41) Ai, are os ] 
ee An,h. Pn ah ( (k+1)¢*1 ) : 


[Ro 
-> (Anghya — Akshy1) — Akshy1) ~E PR-F-A 


ict PR tad ale 
re PI ie Se: n Pn (ener } - 
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Proor: Write gn = y" pn. Now by Abel’s transformation 
n 


0<8 {n,k) = : g (Apsh41 —Akyny1) (Qu-k-1 — Qp-*) 


n 


371 


pay te a Ap (Apshs1) Que — (Anyhgg — Ak+hy1) Gn-k 


n 


ee) (Ap.h+2 = Ap.41) Qu-k — (An,h,2 — Ak4h41) qn_k 


Bak 
n 
< =. (Auinya — Avant) Jp-k 


as gn > O and non-increasing. 


As (R, An, 1) is regular, it follows that 


1 
Ra = PRS > (Ap i42 an Av4h1) p-* 
P= 


converges to a limit asm — co. Now since Pn — ce (by Lemma 8 (i)), 
(4.11) that 


B (n, k) < Rn _ , 
Anihyy Pn Pn att) 


as n — o°, for any fixed k. 


Now by Abel’s transformation, (4.11) and (4.12), we have 


— 1 
Spas >: ne ( oan F; 7, ) k) 


n=kt1 


co 


1 
Aush,2 os Ap+hi qu—-k > An Cer Pn 


(os) 
(Aushs2 cow Apshs1) qu-k 
Apyhyt Py 


co 
< > ( 
bk n= 
bok 


oo 


ag, Sk Fr 


(4.11) 


we have from 


.. (4,12) 


) 
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by Lemma 10 (ii). 
Next, 


h 
had A. Pn_k_1 


1 
sl dS = °| ee) 


n=k+1 
by Lemma (10) (ii). 
This completes the proof of Lemma 12. 


Lemma 13—Let h be any non-negative integer such that vy” pn € DV and let 
An > 0 non-decreasing and 


n (Angi — An) = O (An). 
Then 





chia > ( Akshay Ansh41 ) | Ai e (n, K) 


n>k 


=e (erry 


Proor : Since 


iss Avtt Pn_k p ‘i 
AL 20%) = |) ee eee 
k a) Pn Pn Pn- re ee 


and since by Lemma 8 (iii) 


Pn _ ». 
Pn-1 = 0( n ), 


It follows that 


Sts 





x 


ye (be + hr) 


eg esa 
(k = | etl Ak 


by Lemma 12. 
This completes the proof. 


Lemma 14—Let An > O and non-decreasing. Leth bea non- 


negative integer. 
Then : se! 


l 
n’*l Xn Anti ( =) = 
n tn O (1) 
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implies that 


I 


nr An Ar ( == OL) for r = 0,1, 2, :.. @—-1. 


Proor : This can be proved as in Chow? (Lemma | 3). 


5. Proor or THEOREM 4 


By Lemma 5 (a), it is enough to prove that Jk is bounded. By effecting h + |! 
times successive Abel’s transformation, we obtain : 


n 
8 “4 Q(n, v) cv & 
v=k 


n 
en SS Nahe ( Q (n, v) ) o('+1) 
Ay 


v v-k 
v=k 
hil n \ 
h+ 1 ) SS ( hel ) Ay O(n oft) 
-> ( C4 *y Avr REMI Toe 
r=0 v=k 


We now consider the terms r = 0 tor = h together and the term r = h+1 
separately in the last sigma and obtain 


ato, 








where 
fora) n 
] hat (h41) 
ay’ (ai Oe 
ear > 1 >, Gan he OOM ae | 
n=k v—k 
h 00 = ; 
h+1 >: we : ) 
= > ("Fare D1 D8 ae 
r=0 n=k yak 
(441) 
x Ay A (ny) cy |: 
Writting 


a 
Avehsl1 ¥ Av phy Ansh41 Anghst 
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we obtain 
7g 780) poe 
where 
fee) n 
tanya a sa) (aso ) 
Jo = Xe Pk Sy | Sal Avh41 An4h41 v (n, v) 
n=k v=k 
A est 
fora) oo 1 1 
2 (h+1) ( SS a ) 
O (1) Ak Pk SS Cy-k Av4h41 Anyh+1 
n= aaah 


x Wee Q (n, v) | 


(+) 


—k 
v=k 


by Lemmas 13 and 9. And 


co n 
(12) _ u > n+ (041 
Ji = ARPE >a | Res Q (n, v) one L. 


n=k v=k 


It follows from (4.5) that 
n 


cM) AM! Pay = 1 (n> k) 





vrk 
Vek 
and so 
nr 
(R41) , hod 1 l 
rea aed 9 O(n, = — —- — 
> v-k Pn Pn_i 
Hence 
loa) 
JO?) — 0 (1) Ak Pr by as ( eer 
- Angha1 Pn-1 Pn 
a= 


(equation continued on p. 375) 
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co 


O (1) Pk > ( oe = e ) 


n=k 








= O (1). 
Therefore th is bounded. 


Lastly 


h fora) 
(2) h+1 h Pi 1 
J. vt » ( r ) Ak Pk >: fone A, as ( Aver ) | 





1) v=k 
co 
x1 > 1a ae, 
nv 
ioe] 
| hy1- l 
=0 (1) Pe > hac fale Sal a Na | ) 
( ) — Cy_k (v + 1) ) Avge | 
fore] (h+1) 


u-k 


by hypotheses and Lemmas 11, 14 and 9. 


This completes the proof of Theorem 4. 


6. FURTHER LEMMAS 


We need the following additional lemmas forthe proofs of Theorems 5 and 6. 
Lemma 15—Let h be a non-negative integer and let vy? pn € DH. 


Then 





= eS . 
si ny (n +h) Pn Ol Pp ) 


n=k 
7 (hy2) ; Senin) 
Proor: Since ¢c,, is non-decreasing we obtain : 


2k oo ( A. pnt | Say 


um -(> +> \-ornr 





(equation continued on p. 376) 
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(h42) 2k oo 
2 Ss A» pate fOr 
<asppe D, Me Pe 
(kK+h) Pk mar ee 


h 





eee AREA EE 
(n + h) Pn au Pn-k 


(2,2) 
Ck vy" ? pk x SS (n — k + 1) pn-k 
<a bee (n + hi) Pn Pn-x 


n=2k+1 





K < ae eae 

——— ee a ee <a —— 

< Pk + >) (55 Pn-k Pr 
n=k-1 


by hypotheses, Lemma 7 (vi) and Lemma 8 (v). 
This completes the proof. 


Lemma 16—Let h be a non-negative integer and let (An) be a positive non-decre- 
asing sequence such that : 


nitl Ah*1 An = O (An). 
Then 

n+l Artl An = O (An) 
for 

pie: Oo TZ 4t aels 


Proor : It is enough to prove the result for r = A — 1, since the general result 
can then be obtained by repeated applications of this case. The proof is based on 
one verson of the ‘discrete’ analogue of Taylor’s theorem with remainder and is given 
by the following formula, valid forO < k <n: 


h n_1 
Gee 5 aes ) Aran a > Ca ae ) arta, 
r=0 v=n-k 


(6.1) 


This is easily proved by induction on h = 0; and if in the second term in (6.1), we 
substitute 


n-] 
AMT Ay = AMT An + ZT Abt?) 
b=y 
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we obtain, after simplification, the same formula as (6.1) but with 4 replaced by A + 1. 
This establishes (6.1). 


Note that if c is a constant such that 0 < ¢ <1, then for0 < k & cnandn—k 
<v Qn — l, we have vy > anwherea« > 1 — c > 0. Hence 


Ay A A 

h+1 = Ot die Oa ee : ee 2 

AMT Ay o( 5 )=0(—h-)-0( ae ) 

as (Av) is non-decreasing. Since the sum of the coefficients in the second sum in (6.1) is 


k+h 
( ee ) = O (kht+1) = O (n**1) 
it follows that, uniformly in 0 < k < cn, we have 


h 


Pit iis , (> Es ) Report 6.2) 


r=0 


Now let d be a constant with O < d < ; , and let k =[nd], It follows from (6.2) 
that 


h 
> ey (fF )1e = 0.09 63) 
v=0 


since the coefficients are constants. Now the sum in the left of (6.3) is the Ath diffe- 
rence of / (n, vk) (regarded asa function of v) taken for v = 0. The terms withr < h 
in the sum defining / (n, vk) are polynomials in v of degree less than h; their Ath dif- 
ferences are therefore 0. The terms with r = hisa polynomial of degree Ah, the 
coefficient of v” being 


Pa ANN, 
(—-1P hi 


its Ath diflerence is therefore 
kh AM An. 
Hence (6.3) reduces to 
An An 
nan = 0 (a2) = 0 (4) 


by definition of k. Hence the result. 





Lemma 17—Let An > 0 and non-decreasing. Then the hypotheses (3.4) and 
(3.5) are equivalent; i.e. 


nel An Arti ( =) = 0 (li) 
An 
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if and only if 
n+l Ah*l An = O (An). 
Proor : Suppose that (3.5) holds. Then by Lemma 16 
n+l Art] Xn = O (An) 


for 
pom L 2 0 2 ees A 


Write he for any product of terms each of the form An,» where bisa constant in- 


teger. The values of b for different factors of the product may not be thesame. We 
use a similar notation for (A’ An)(®). When a = 0, this is taken as meaning 1. 


Since 
1 A An 
. ( i) age ph 
we can verify by induction on hthat A?*1 ( aa ) is the sum ofa finite number of 
n 
terms of the form : 
(AAn)(%) (A2 An)(@o) .. (A*+1 An) @n41) 
(0) 
An 
where aj, ao, ..., are non-negative integers such that 
a, + 2a2 + 3a3 4+... + (Ah + ll) any =At+ 1 
a4 + ao+...+an4,+1=b, 


It follows easily that (3.4) holds. 


The converse implication may be proved in a similar way using the formula just 
stated with An replaced by 1/An. ‘ 


7. PROOF OF THEOREM 5 
By Lemma 5 (b), it is enough to prove that Jz is bounded. Now 
he< id 4 1@) 
where 


Kj n 


co 

7) pt ] 

k . an Pau _ Pp > Cv-p Pn_v Ay | 
“= 


B=0 v=k 


MODIFIED MEANS 379 


n 


oo 
] 
he = Pk-1 eae >, Ay Pn-y Cy-k 
n=k 


v=k 


By Abel’s transformation A times. 


n 
xX Ay Pn-y Cy-u 


vk 
= Sa. Pn_y) Pade —— Paha 1 Ak Pn-k 
v=k 
n h 
= > Ne (Av pn_v) fee = > c ae AG: (Ak pn-t). 
vok r=] 
Hence 
(1) (11) (12) 
I, </ k + Bs 
where 


00 
6 fa hor 
goed (15% Sake Sern 


r=0 p=0 n=k 


h 
(Ar pn-v) ou | 


h k-1 20 oo A’ pn_v 
hr v 
2 NDS 2S Gate Sera 
r=0 P=0 vek oe 
h hi co [o | | (So Ager | 
= O (1) Ss ( r ) Sy Pp Ss) yr Av 
r=0 P=0 vek 
k-1 0 Ic oe 





\ 
ie) 
M 
> 
oe 
“4 
+ 
J 


b=0) vak 
k_ 
1 : (”) l 
= O(1) DS Pp Crag 
p=0 bert 


(equation continued on p. 380) 
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k1 
l (241) 
=i O41) => Pp CHS 
#=0 
1 
= O(1). —7,— kh = 0 (1) 
by Lemma 8 (iii), hypothesis (3.5), Lemma 7 (iv) and by the fact that 
n 
Dd pare? aa 
k=0 


which, by virtue of (1.21), follows from the identity : 


(E pn xn) (Sc) xn) = Og? xn (7.2) 
n n 


Before we consider hai » note that by hypothesis (3.5) and Lemma 16, we have 


n (An a An41) =O (An) 


and this implies that 








An ~ Xn41. 
Hence 
eet | (7.2) 
as n — oo, for fixed @. 
Now 
k-1 h 
i) e S; rn = | BS Pp > cee AT" (Ax pas) | 
n=k p=0 r=] 
oo , k-] ee | 
= S Acta > Ss) es ) rig 
n=k w= r=1 6-0 


x ( a, * eee) A pn-p | 


r-1 


h co 
r— | ] Rae 
<a 9 i >s4= (A, * Akye) 
§=0 n=k 


(equation continued on p. 381) 
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. A; pn_ Sh Coe 


F=0 


co 
eee. l io 
ow S ( 6 ) An Pn | A, Akye | 
n>k 


r=1 06=0 





x (Ay Pn_k) kt-1 


rl ar 


ow> Ceaare: Aerie doko S 


at Pah nk 


ll 
M 
M 


| r-l1 
r ee Kinde. Nig Arse | 


ee 
0 Ak 


h f-1 


on SS (151) B20 


r=1 e=0 


by Lemma 10 (ii), (7.2), the hypotheses, and the fact that for r > 0, 
5 pk tae = Anis F 
k=0 

(This identity follows from an identity similar to (7.1)). 

Thus 


1 < 18) + 12” = 0 (1). 


Now we consider i?) . By making A + | times Abel’s transformation we have 
nN 
x Ay Pn-v Cv_-k 
yrk 


h n 


3 ( = S (Ay *** Aver) (A, Pn) Se ) 
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Hence 





h fore) n 
(2) a h 1 h_r41 
= > (.) Dare | Oe ae 


r=0 n=k vok 


h 
x (AS par) ct? | 


>t ) Pe S 1 at ae PS — 





A Sabi Aver | Pr 
-on > | +) > 
v=k 
h co o M+) 
h uk 
-00 S(t) Sie -00 


r=0 v=k 
by Lemma 10 (iii), Lemma 9 (a), Lemma 16 and hypotheses. 


We note that the result of Lemma 9 (a) is valid with the assumption (3.6) in the 
case h = 0, and this is not needed in the case /) > 0. 


(2) . 


Now we consider /," in the case h = 0, under hypothesis (3.7) in place of (3.6). 


Now 
n 
0 (n, k) — a Ay Dn-yv Cy-k 
: ° 
= x A Ay = _ Pa- Bb Ch-k + Anyi Sn,k 
where 


n 


Snk = SS Pn-p Cu-k = { 


Pook 


l(n = k) 
O(n # k), 


Now by hypothesis (3.5) and (3.7) (in the case h = 0) we obtain, 


sv=0(%)-o( +) 
v n 
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Hence by Lemma 7 (vii) 


n 
6 (n, k) = O(1) pax S, HAAG coro tclets tok 
v=k 


r 
= O(1) <* pax ASE As 


n-k 
Hence 
(2) 
W Dak Co. 
(2) n-k Pr. dk 
k (1) Pe > n Pn + At. Pk 
n7k 
= 0 (1) 
by Lemma 15. 


8. PROOF OF THEOREM 6 


This follows by combining Theorem 4 and Theorem 5, taking note of Lemma 17. 


9. COROLLARJES 


We know that Cesaro mean and the modified Cesaro mean with weight m are 
equal. In the following corollaries we examine special sequences A for which Nérlund 
and modified Nérlund means are equivalent. 


From Theorem 4, we obtain the following. 
Corollary 1—Let vy" pn © AZ and let 
An = n® (log n)*. 
Then 
|N’, p, Al | NPI 
in the following cases : 
(i) h > Oand either B > 0, Srealor@ = 0, 3> 0. 
(ii) 4=0, B> 0, 8 real 
In particular 
(C’,a,A)>|C,a| 


in the following cases : 


I 
— 
o7 
V 

Co 


(i) «> land either B > 0, 6 real or B 


(ii) 0<a<1,8>0,5 real. 
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Similarly from Theorem 5, we obtain the following. 
Corollary 2—Let 9" pn € GD and let 
An = n® (long 7)§. 
Then 
|N, P| >I1N, pAl 
in the following cases : 
(i) h>0O, B> 0,5 real 
(ii) A = Oand either B > 1, 6 real or (3.6) holds. In particular 
1 ep eC eA f 
if8 > 0, d real. And 


1 1 
| Pig kl o eeare a 





yA 


if B > 1, 5 real. 
Combining Corollaries | and 2, we obtain the following. 
Corollary 3—Let ¥” pn € ADT and let 
An = n® (log n)8. 
Then 
|N’,P,A|~ 1, P| 
in the following cases : 
(i) h>O, B > 0,5 real 
(ii) A = 0 andeither B > 1, 6 real or (3.6) holds. In particular 
era. Fer ha) ee . 


if 8 > 0, 5 real and 


Noe ee ON eee 
n 


l 
nk + 1 
If 8 > 1,8 real. 
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SCATTERING OF A COMPRESSIONAL WAVE AT THE CORNER OF A 
QUARTER SPACE 


NARINDER MOHAN AND P, S. DESHWAL 
Department of Mathematics, Maharshi Dayanand University, Rohtak 124001 


(Received 6 May 1988) 


The problem of scattering of a compressional wave at the corner of an 
elastic quarter space (x A 0, z A 9) has been discussed inthe present paper. 
One face of the quarter space is free and other is assumed to be rigid not 
permitting the displacements across it. The technique is due to Wiener and 
Hopf. The scattered field possesses the character of transverse cylindrical 
waves. The numerical computation for the amplitude of the scattered field 
exhibits a sharp fall yvcrsus the small values of the wave number. 


INTRODUCTION 


A problem of special interest in seismology is the problem of scattering of an 
elastic wave at a corner of an elastic medium as it leads toa phenomena of a surface 
wave entering into another medium after travelling through a medium. Several 
authors* 611 have studied the problem of Rayleigh wave propagation in elastic 
wedges both theoretically and experimentally. Recently, Momoi’ has discussed the 
problem of scattering of a Rayleigh wave in an elastic quarter space. He has extended 
his study to Rayleigh wave scattering due to a rectangular mountain’. The three- 
dimensional problem of scattering of waves in an elastic quarter space has been con- 
sidered recently by Gautesen®. Both Momoi?’8 and Gautesen? have studied the pro- 
blem of scattering of surface waves in an elastic quarter space using rigorous 
numerical computations to get approximate results. Deshwal and Mann! have at- 
tempted the problem of scattering of Rayleigh waves at the corner of an elastic 
quarter Space to obtain exact results using the technique of Wiener and Hopf. 


In the present paper, it is being contemplated to study the problem of scattering 
of a compressional wave at the corner of a quarter space. Most of the mountains are 
deep-rooted with their bases in the solid mantle of the earth. They are assumed to be 
rigid forming the rigid boundary of the present problem. They peep out of the earth 
and scatter the seismic waves. Physically, mountains within the earth form the rigid 
boundary and the surface of the earth is the free surface of the quarter space. The 


Fourier transformation and the function theoretic technique due to Wiener and Hopf 
is the method of solution. 
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Basic EQUATIONS 


The problem is two-dimensional. The waves propagate in the zx-plane. The 
x-axis is in the free surface and the z-axis being along the rigid boundary with origin 
at the corner of the quarter space x >0,z > 0. The medium is homogeneous, 
isotropic and slightly dissipative. A time-harmonic two-dimensional compressional 
wave is incident at the corner (Fig. 1) and gives rise to the reflected and scattered 
waves. 





Fig. i. A solid quarter space with a rigid boundary. 
The incident wave is | 
di (x, z) = exp [-—ik (x cos 6 + z sin 6)}. ak 


The total potentials in the medium are 


bt (x, Z) = gi (x, z) + $ (x, 2) 2} 

wt (x, Zz) = (Xx, Z). (3) 
The wave equations are 

(v2 + k2)¢ =O = (v2 +k) ¥. ...(4) 


k = ky + ike, kk’ =k, + ik’, are complex wave numbers with positive imaginary 
parts. The displacements (u, w) are given by 


Bare ON Ob. BU (5) 


SOx. éz ’ dz ox 


BOUNDARY CONDITIONS 


The conditions on the boundaries and at distant points of the quarter space 
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are 
(i)u=O0=w,x=0,z7z20 ...(6) 
(ii) Pez = 0= Pe, Z = 0, 5 0 PAYS, 


where pzz, Pzx are the normal and the shear stresses. The conditions (6) imply that 
¢: and $+ satisfy the Cauchy-Riemann equations and are harmonic functions. 
Equations (4) will result in 


dt (x, z) = 0 = Ue (x, z), x = 0,z 2 O. | ~«e8) 


The half-range Fourier transform 
a co 
$+ (p,z) = J d(x, z) et?* dx, p = a + if ...(9) 
0 
is analytic along with its derivatives in the region 8 > — dof the complex p-plane if 
for given z 


| d(x, z)| ~ M exp (—d|x|),M,d> 0. me, 


The Fourier transform of v (x, z) and its derivatives with respect to z are analytic in 
the same region. 


DISCUSSION OF THE PROBLEM 


Let us take the Fourier transform of the first of eqns. (4) to obtain 
d2 os b¢ 
ea, Bad Patel *.. 5 
(5 a Yeu (p, z) ( ax ), ip ()o -«(11) 


where Y = + (p* — k?)1/2. The sign before the radical is such that the real part of 


Y is always positive for all p. The subscript 0 denotes the value at x = 0. The first of 
the conditions (8) simplifies to 


(¢)0 = — (¢i)o = — exp (—ik zsin@). ~ mti2) 


A complete solution of (11), which holds when z > co, is 


rE —_ Sig 2ip exp (—ik z sin 6) 
= = = og Uheee 
$+ (P,z) — $4 (—p, 2) = A(pe oh ae atlas 


Putting z = 0 in (13) and in its derivative with respect to z and eliminating A (p) 
between the resulting equations to obtain 


$.(p) — $4(-p) = — | +, (p) — $', (—p) 


+ 


2pk sin 0 ] 2ip (14) 


p2-k2 cos? 6|  p2—k2 cos29 ° 
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The notations ¢. (p), ¢, (p) are used for ¢, (p, 0), ¢, (p, 0) etc. Equation (14) 
is a Wiener — Hopf type functional equation to be solved for two unknown functions 
Ps (p) and ¢, (—P). 


SOLUTION OF THE WIEBNER-Hopr EQUATION 


Equation (14) can be written as 


et! Gyn ak i 
= ke 80 8s 
p+ (P) + Y ae Y(p+k cos 6) ra $4 (—P) 


ays F (pak cos 8) 


(15) 


The left-hand member is analytic in the region 8 > — ks cos @ and the right hand 
member in B < kg cos @. By analytic continuation, they represent an entire function. 
The points p = + kK are excluded by the branch cuts. Each member tends to zero as 
|p| > oo. By an extension of Liouville theorem each member is identically zero. 
Hence 


,, (p) i : 
a iY + k sin @ 
¢,(p) = ~ a ae VUTEC ea ae .. (16) 
Similarly 
v. (p= —¥, (PBS = + (P? - k'2)112 (17) 


the choice of signs for 5 is same as foter.. 


faking the Fourier transforms of the conditions (7) and using (16) and (17), it 
is obtained that 


(2p2 — k’?) v4 (P) 2ipk sin 0 
P P+ 5 p-k cos @ 
BRIE Poesine tay 
oS 3 saad rip ¥, (») = [ Y(p + kcos @) 
fF (2p — k’). 9 
+ pee | OP St: 


These equations are solved for V’, (p). p is ‘changed to —p to find the value of 


. (p) — ¥, (-p). Then 
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i mete ke (py = Wh (—p) |exp (—8z) 


_ 8puk cos 6 (2p? — k’?)(¥ — ik sin @) ob: 
y: (p2 —k2 cos? @) F (p) ‘ 


...(20) 
where 
F(p) = »[(2p2 — k’'2)2 — 4p2 Y8}. Het Ob 
VARIOUS WAVES 
The potential function (x, z) is given by the inverse Fourier transform 

ie a) +iB 

bonz)— 9 | Ty (pz) et dp, x > 0 (22) 
7 
— CO+iB 


where —d < 8 < d. The factor exp (—ipx) = exp (ixx) exp (3x) in (22) vanishes 
as p — — co in the lower part of the complex plane if x > 0. If the contour of inte- 
gration is chosen in the lower half of the complex plane, where /, (—p, z) is analytic, 
then 
oc +88 
bona = 5> | fhe (2) - Cr, dete ap, ...23) 


—oo +iB 


The value of expression within the brackets is obtained in (20). (23) is evaluated along 
a closed contour in the lower half of the complex plane with p = -—k, —k’ as the 
branch points, indentations around p = +k cos @ and satisfying the conditions 
Re (Y) = 0 = Re (5) (Fig. 2). The conditions according to Ewing et al.2 imply hyper- 
bolic paths for the branch cuts at p = —k, —k’. Indentations around P= +kcos @ 
contribute the reflected transverse waves 


D (k) exp ( -ikx cos @ — 8’ z) and —D (k) exp (ikx cos 9 — 5’ z) ...(24) 
where " 


D (k) = 4 » k* sin @ cos @ (2k? cos? @ — k'2)/F (k cos @) +0a(25) 


and 5’ = (k2 cos? 9 — k’2)tl2, These waves cancel each other on the rigid boundary. 
The compressional waves reflected from the boundaries are 


— 4 exp [~—ik (x cos 6 z sin @)] and — 3 exp [ik (x cos @ — z sin 6)]. 
++.(26) 


These waves together with the incident wave vanish on the rigid boundary satisfying 
the boundary conditions (8). 


The scattered waves are the contributions of the integrals along the branch cuts. 
Along the branch cut at p = —k, real part of Y is zero and imaginary part of Y has 
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Fic. 2. Contour of integration in the complex P-plane. 


opposite signs along the two sides of the branch cut. 5 is unchanged along the cut. 


The integrals cancel each other along the branch cut atp = —k. Thus the scattered 
waves are not of the compressional type. For the contribution along the branch cut 
at p = —Kk’, let us take p = —k’ —iu, u being small, then 

§2 = (—k’ — iu)? — Re = 2k, ik, ) -—uv 122) 
Along the cut, real (5) = 0 and 32 is real and negative. Thus 5 = +i0', 
8’ = (2k, 3 Poa id 0 and u varies from 0 to °° along the cut. Evaluating (23) 


along the sides of the branch cut to find 


es ep 





ie 2 ore - 
be 2) = —p— || fo 2) — % CP ete 
0 


= 10, (Pp, 2)'— pb. (—P; Ze’ | e~¥® du _..(28) 
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{( (u) sin (4) 24 uz + N (u) Vu cos (af 2. u z)] 
0 


(ka + ud) x 4, (29) 








where 


8pk cos 6 (Y — ik sin 6) p2 (2p? — k’?)8 .. (30) 
Me = H(p) (p? — k? cos? 8) 


3283 Yk cos 6 (Y — ik sin @) (2k 2)1!? u2 (2p? — k’) 


wel) 
n H (p) (p2 — k? cos? @ 


N(u) = — 





H (p) = p? [(2p? — k'?)4 + 16p4 Y? 8°]. el 3 


Since u is small, M (0) and N (0) are retained. Following Laplace integrals are used 
for evaluation of (29) Oberhettinger!® : 


co ———— af 28 - k, zs 
[sin (,/ 2k, uz | e- ut dy = exp (- ac ) pea (2 
0 


~ 9x82 12 








oo Vm (x—k; z2) ke 2% 
| Ju cos ( 4] 2k; uz ) ewe du = ee ee exp (- oe ) 
0 





...(34) 
Further, when x > > z, then 
r= (x? + 22) 1/2 = x (1 + 22/x2)ll2 = x + 722/2x, ...(35) 
From (39), it is obtained that 
=e «| 2k sin 8 
ta (x 2) = is craat (0) + say 
cos 8 — k’ r sin? p 
+ N (0) ——____,______ ...(36) 


r cos 5/2 8 J 
where x = rcos 8, z = rsin $8. Thus the scattered waves in (36) are transverse waves 
behaving as cylindrical waves. Their amplitude is of the form exp (—k, r)lV/r which 


decays exponentially as the distance r from the scatterer at (0,0) increases. In the free 
surface (z=0), it can be seen that (33) and (24) behave as O (x~3/2), Close to the 
scatterer as r>0, the amplitude of the scattered wave behaves as O (r°3!2), 
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It is interesting to note that the scattered waves are transverse and not com- 
pressional. For the far-field when x >> z, the amplitude of the scattered wave 


behaves as exp (—k, r)/V a r = (x? + z2)1/2,r being the distance from the scat- 


terer at (0,0). The wave is a cylindrical wave. It dies out exponentially as it moves 
away from the scatterer. For the near-field as r>0, the scattered wave is of the form 
O (r-3/2). It is dominant near the scatterer. The amplitude of the scattered wave is 
plotted versus the wave number (fig. 3). For Poisson’s solids, when « = +/3 8, @=2/3, 
the amplitude in the free surface (z=0) of the quarter-space falls off sharply as the 
product of the wave number and the distance increases slowly. The amplitud has the 
value 36.60 when kz x = .1 and falls to the value .025 when kp x = 1.8. There is no 
scattered wave when @ = 7/2, i.e., when the wave is incident on the free surface 
parallel to the rigid boundary. The waves in (24-25) (both transverse and compres- 
sional) reflected from the boundaries of the medium satisfy the boundary conditions. 
The scattered waves together the waves contributed due to the branch cut at p = k’ 
also satisfy the boundary conditions. 


1.5 
1-4 
1.3 
1-2 


7 0 6 
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Fic. 3. Amplitude of the scattered wave. 
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The residues at the poles p = +/pn contribute 


8 pn k cos 6 (2p" — k’*) (Yn — tk sin @) 
+ oii e Fn? etiDn® 
(p; — k? cos? 6) F’ (+pn) 


where Yn = (p. — k@)ll2, §n = (p. ~ k’2)1l2. These are the Rayleigh waves as SV- 


type waves and P-type waves. These are surface waves. Their amplitude does not 


depend upon the wave number and is a constant multiple of the amplitude of the 
incident wave. 
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Green and Lindsay! have given the linearized temperature-rate dependent 
thermoelasticity theory. Using this theory the problem of thermoelastic 
longitudinal vibrations of an infinite circular cylinder has been solved and 
the results have been obtained in terms of potential functions. It is interest- 
ing to note that due to the temperature rate dependent theory the amplitude 
of both the elastic and thermal waves are higher than that of conventional 
theory. We further observe that if we put y = y* = 0 inthe results of this 
paper we arive at the results of the conventional theory of thermoelasticity. 


1. INTRODUCTION 


The thermoelasticity theory which includes the temperature-rate among the 
constitutive variables has aroused considerable interest in recent years. This theory 
which was developed by Green and Lindsay! and by Suhubi? is a generalisation of 
the conventional coupled thermoelasticity theory? and predicts a finite speed for the 
propagation of thermal signals. Several problems revealing interesting phenomena 
characterizing this theory have been considered earliert®. Because of the experi- 
mental evidence available in favour of the finiteness of heat propagation speed’, these 
studies are of practical relevance too. 


The problem of thermoelastic vibrations of a circular cylinder was discussed by 
Chadwick8 and the solutions have been obtained in terms of potential functions. In 
this paper an attempt has been made to study the problem of thermoelastic longitu- 
dinal vibration of an infinite circular cylinder in the context of linearised tempera- 
ture—rate dependent thermoelasticity theory and the solutions have been obtained in 
terms of the potential functions. It is interesting to note that due to the temperature 
rate dependent theory the amplitude of both the elastic and thermal waves are higher 
than that of conventional theory. We further observe that if we put y = y* = 0 in 
the results of this paper we arrive at the results of the conventional theory of ther- 
moelasticity discussed by Chadwick®. 


2. Basic EQUATIONS 


In the context of the linearised temperature—rate dependent thermoelastic 
theory of Green and Lindsay!, the equations governing the displacement vector 0, 
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the stress tensor o1j and the temperature deviation @, in a homogeneous and isotropic 
solid, in tensor notation are given as follows. 


otj = A Sag uke + p (Utg + 49,4) — 7 (@t+y 8) 847 zt 2e1) 
pPO=pyO+ (Act p) grad div 0 — 77 brad (9 + ¥6) ate) 
k 92.9 = Pee (8 + y* 6) + a i,k. (2.3) 


Here P is the mass density and A and »p are the isothermal Lame’s constants given by 


a ee (2.4 
eee RS Ey 2 (1 + vr) mie 


in which Er and vr are respectively isothermal Young’s modulus and isothermal 
Poisson’s ratio and 


i Fase) (KS) 


is the coefficient of volume expansion. 


The isothermal compressibility is given by 


ct 0 uk,k es 1 — 2vr 
ef a 3( Oo okk ), ™ 3( Er ) ay) 


from which we see that Er, vr as their suffixes indicate, are the isothermal elastic 
constants. Further Ce is the specific heat at constant strain, k is the thermal conduc- 
tivity, 4 is the initial uniform temperature and yY, y* are thermal constants charac- 
teristic of the theory. A superposed dot denotes partial differentiation with respect to 
time ¢. It is assumed that the body forces and heat sources are absent. Expressing the 
displacement vector U as the sum of irrotational and solenoidal components 





C= yv¢+curlA. (2-7) 
Equations (2,2) (2.3) are found to be equivalent to the set 


024 Er (1 — yr) Er 


eg re ee 5 
or? (1 + vr) (1 — 2vr) vi ¢ 3P (1 — 2vr) (06 + y @) 


2 pe Er 8 « 0 
ky26@ = PC, * ieee See ee 
VO = PCO + 8) + ao veg (2.9) 
ae ti pO 
or2,—~=<“<i—i‘«‘SW (1 + vr) V~ 4. lea tun 


From equation (2.8) we observe that the velocity of longitudinal elastic waves in a 
medium with zero coefficient of expansion is given by 
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a Er (\ —_ vr) 1/2 
Vr | wey (Le ave) : Rate. it) 


We shall refer to Vr as the isothermal velocity. Further the velocity of transverse 
elastic waves? is given by 


= 1/2 
aa e p(l + vr) ] (2.12) 


3. THERMOELASTIC VIBRATIONS OF A CIRCULAR CYLINDER 


We denote by (r, 9, z) cylindrical polar coordinates referred to the axis of the 
cylinder and make use of potential functions. In longitudinal disturbances the tangen- 
tial component of displacement vanishes identically and since, in addition, field 
quantities do not depend upon the angular coordinates no confusion will arise from 
the use of @ to denote the temperature perturbation. The vector potential A is of the 
form (0, ¥, 0) and the displacement components are given in terms of the scalar 
functions ¢ (r, z, t), v (r, z, t) by 

op oy se 


ur = oar so Oz” us —— 0, he 7s eof (r y). pa(3.)) 


The thermoelastic equations (2.8) — (2.10) now take the form 














a2 ¢ ee 1 06 @2°§ a. tian 
ae ee = eT chore te Da Py ee 76) ...(3.2) 

a0 8 ( &@¢ 1 8d aes) 
pC. (9 + y* 6) + rz | a ee rem 

62 6 1 06 a2 @ 

nc ( pratt riira aos ) ---G.3) 
al AUER ab Aaa 4 ) 3.4) 
or aa “Or2 a r\ or r2 ¥ 6z2 ( 

Inserting into equations (3.2) — (3.4) the solution 

(0, 6, ¥) = (0 (r), o Cr), ¥ (r)) Orn. rsh 3s9) 


In these expressions y and © are, in general, complex quantities so that the wave 
length is 27/Re » and the period 27/Re ». 


A . Pat n 
C2) Re EC ...(3.6) 
dr2 r dr y? pxr y2 


d29 1 dé Ee GQ —iytey—w]e 
at r dr ~ ; 


(equation continued on p. 398) 
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imaGo d2 6 1 de 7 
(Sf+2ge-78) a7 





kh 
ei yt el 114 
Se +[( - - #)- 4] 4-0. (3.8) 


Solving these equations, we get 
b= Ado (raf & — 1) + Blo (ral & — 2) exp ali(nz—o1)) ..3.9) 
PXr ve 


ene rates a 6s pas ei ) Jo( raf & — 7 ) 
> 


ra lta Le ) Jo ( 4] ——— n?) | exp (i(nz — wt) 





polonu 
pe Ch (ral 2s — 9?) exp (i (nz — wf)) honL as 
where e ; z are the roots of 
v2 
iw3 °C. 
+ -———. (1 — ioy*) = 0, (3.12) 
= Ri & S 
kV". 
The equation corresponding to (3.12) from Chadwick® is 
2 ? 
ga — | = + ac (1 + )| cle 
ye 
iw? PC, 
a oenraee ee. 
= ...(3.12A) 
KV. 
: a* do 
and 
w2 
Gm 
3 3 .-.(3.14) 
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The Bessel function Yo, ¥1 of the second kind are excluded from the solution by the 
requirement that wy, uz, 6 shall be finite at the axis. 


We find from (3.12) and (3.12A) that due to the introduction of the temperature 
rate theory the values of EY + 3s and c & are reduced by amounts of iw (ey + y*) 


and iwy* respectively. This implies a reduction in the values of + and As which by 


equations (3.9) and (3.10) cause increase in the amplitude of both elastic and thermal 
waves. 


The constants A, B, C are found by substituting (3.9) — (3.11) into boundary 
conditions which, since the surface of the cylinder is free from mechanical and 
thermal constraints, take the form 


Orr = 0, yz = Oat r =a. aoa ia EY 


Equations (2.1) and (3.15) give 


Dae et Aad 1 8¢ 824 o(1 a 1 
orn. = ev’ ( ore + —_ SOs az2 — 2pV, ( or 











| 

| 

eS) = 54 v0) = | 

+ 622 am erez tr (6 2 | 6) kee 0 | 

Pee ee a =a 0. f 

sla JA ( ac; 6z 35 er2 r Fs Or r ez2 : | 

Thermal boundary condition is | 
m [atr =a 

6 
ér jy ...(3.16) 


where hi is the heat transfer coefficient of the boundary r = 4. 


The three homogeneous linear equations connecting A, B, C are found to be 


[0-34 








2 
= Jo (anHi) + = Ji (a4 mh) | 


salt ~ 7) 


2 
gei0 








2) jo(an He) + 2 (an Ha) | 





I — 
=: 2iCHs | Jo (a n Hs) > anHs3 Ji (ants) | = 0 SACRED 





a Beg ( 
2i [AH1 Ji (anf) + Be Ji (anH2) — © (2 pas ) Jy (a0H3) 
& 


...(3.18) 
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A Jo (anHy) — WM Jy (ant) | ( pa os Hy )] 





2 
Ne 12 
A 2 

+ a = Jo (anH2) — He Jy (anH2) | ( ele A, )=0 

) V > 42 

r 4 
fe cts 

where 
Sec 

H; ati ae joe? 1, 253; Bes 9-1) 


Eliminating A, B, C between equations (3.17) — (3.19) we obtain 


(2 - 


@2 








2 
P ) Ji (ands) H( w73 +1-— A ) Jo (ay He) Ji (an Ay) 
an "2 Vi 1 


@2 





— Ho (#2 +1— ) Jo (an Ay) J, (anHe2) ] 


2 
ae 72 


2 2 2 
+ 4, Hi Hs (# ital \n (an Hi) Ji (nH) 





2 
x {20H Jo (anH3) — — J; (anHs3) } 
vy? 72 


& 





I 
3|> 
rH 
— 

nN 

| 

N 
i) 
rn 


@2 2 2 2 
(i - ) Jo (arth) 


x Jo (anH2) Jy (anH3) — = 2 anH3J (anH3) 











_ _w2 2 2 
A (onli) x (HY +1- — ) 
y? "2 V" 42 
8 r 2 
x Jo (anHe) Jy (any) — Ho ( #7 +1] — < ) 
y? 7? 
a 
x Jo (any) Jy (anHe) a me GPP 


Equations (3.12), (3.14), (3.20) and (3.21) together determine the wave number the 
wave number 7% as a function of the frequency «. 
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We observe that the form of equation (3.21) is the same as the frequency equa- 
tion obtained by Chadwick’. The difference lies only in the values of Hi and He 


which in the present case depend on y and y* also and if we put y = y* = 0 the two 
equations become identical and the results tally. 


If we replace Vr by Vp and put 


yp2 
Hy = > area aages Fe Hy = co 


=(3:22) 
2 
ae 


where 


v= — 3.23) 
x] 


Equation (3.21) reduces to 


TT 


2 \2 2 \ oe 2 Rie 
G- FY a (nd) a(n) 
4 a Vi 

3 ol V2 
+(2J— -1) Ji (an 4|-—- 1) 

an V 12 p2 

Dp D 


SS 


2 eg 
x { (200 4 = = 1) Jo (an) - 1) 
V V- 


8 


Ebest 
eden ga (a omen 1) |- 0. ...3.24) 
2 y° 


Equation (3.24) is the final result of the Pochammer Chree analysis of the longi- 
tudinal modes of Vibrations of a circular cylinder. 


To find the order of magnitude of the discrepancy between equation (3.24) and 
the exact frequency relation (3.21) we expand the roots of equation (3.12) in powers 
of the reduced frequency %1. From equation (3.20) we obtain 





HN = ‘Lore — 1 [Il + O(%)] et § Ber. Y | 
y2 
Dp 
Pee tt 4) (1+ 0 (2)] (3.26) 
OT) 
- BER Beg 
using the relation Vp = Vr V1 + € 
eCeVr 


© 
‘in? 





@ 
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It follows that 


y2 
Jo (ax) = Jo Rasp -1) + O(%) asXt—>0 
. V 


p : 
Jo (anHz) = (3X)1” learicee} 


V ; Lae 
x exp {Fem +e) (1 — i) +" i \, ... (3.28) 


with similar results for the first order Bessel functions J; (ayH 1) J; (ayHe2). Entering 
these expressions into equation (3.21) it is found that each sider educes to multiple of 
L + O(%) where L is the left hand side of equation (3.24). Thus in the limit % > 0. V 
becomes indererdent of % and h and we recover the classical frequency relation (3.24). 
Since this equation has been shown to differ from the thermoelastic frequency relation 
* (3.21) by terms of order %, the error sincurred by using the Pochammer—Chree analysis 
in the interpretation of experimental results are generally quite insignificant. 


REFERENCES 
1. A.E. Green and K.A. Lindsay, J. Elasticity 2 (1972), 1. 
E. S. Suhubi, Thermoelastic solids. In Continuum Physics (ed. A.C. Eringen) Vol. II, Part I, 
Chap. II. Academic Press, New York, 1975. 
M.A. Biot, J. Appl. Phys. 27 (1956), 240. 
A. E. Green Mathematica 19 (1972), 68. 
V.K. Agarwal, Acta. Mech. 31(1979), 185. 
D. S. Chandrashekaraiah, Proc. Indian Acad. Soc. (Math. Sci.), 89 (1980), 43. 
C. C. Ackerman and R. A. Guyer, Ann. Phys. 50 (1968), 128. 


P. Chadwick, Thermoelasticity— The Dynamical theory. In Progress in Solid Mechanics, 
Vol. I. Amestardam. 1960. 


3 


PNA AY 


Indian J, pure appl. Math., 20 (4): 403-422, April 1989 


LINEAR STABILITY AND THE RESONANCE FOR THE TRIANGULAR 
LIBRATION POINTS FOR THE DOUBLY PHOTOGRAVITATIONAL 
ELLIPTIC RESTRICTED PROBLEM OF THREE BODIES 


V. KUMAR 
M. S. College, Bhagalpur 


AND 


R. K. CHOUDHRY 
University Department of Mathematics, Bhagalpur University, Bhagalpur 812007 


(Received 12 February 1988; after revision 8 August 1988) 


In the present paper we have found the range of values of p and e for the 
linear stability of the triangular points for the doubly photo-gravitational 
elliptic restricted problem of three bodies and it has been shown that 
some resonances of the third and the fourth order exist which will need special 
investigation for the determination of complete stability of.the lihration 
points under our reference. 


1. INTRODUCTION 


Simmons et al.14, studied the existence and linear stability of the libration points 
for the photo-gravitational circular restricted problem of three bodies. They found the 
ranges of the linear stability as well for all the various libration points. Manju and 
Choudhry? studied the non-linear stability for resonance as well as for non-resonance 
cases for the isosceles triangular libration points when only one of the gravitating 
bodies was taken to be radiating and the reduction factor very small equal to that for 
the sum. The authors® studied the non-linear stability of the triangular libration 
points for non-resonance case for the doubly photo-gravitational restricted problem of 
three bodies where both of the gravitating bodies were assumed to be radiating and 
similar to the assumption of Simmons and others in the mentioned paper we took the 
reduction factors to vary from —o< to 1. In another paper® we have studied the non- 
linear stability of the same libration points in the presence of the third and the fourth 
order resonances. 

r we have studied the linear stability and its range for the 
ity of the elliptic orbit. We have also shown 


h order will exist under our range of linear 
bility we shall have to take into considera- 


In the present pape 
various admissible values of the eccentric 
that resonances of the third and the fourt 
stability. So for the study of non-linear sta 
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tion the resonance as well as non-resonance cases. Our prescnt paper may be claimed 
to be a generalization over our mentioned papers in the sense that we consider here 
elliptic restricted problem and over Markeev’s work! in the sense that we have here 
assumed both of the gravitating bodies to be radiating as well as where Markeev has 
taken them only gravitating. We have not investigated the existence of the various 
libration points in detail. However, the existence of three collinear and two triangular 
libration points can be shown. The present paper studies the linear stability of the 


triangular libration points only. A similar study for the collinear libration points will 
be interesting. 


The present problem which was first taken up by Radziavskiil® was studied in 
detail by Simmons ef a/.14, They remained restricted to the circular case and contented 
themselves with the linear stability only. Realising its far-reaching consequences we 
took up the investigation of the problem in full detail. The results are of immediate 
importance in stellar dynamics, but it is not less important for the solar system. 
Regarding the two reduction factors to be arbitrary parameters, the problem reduces 
to one where the gravitational forces may be taken to vary from the classical case 
(x = B = 1) to the trivial case (« = 8 = 0), Negative values for « and 8 can be 
available only for stellar system. Thus we find that the results are equally applicable 
to the solar system as well as to the stellar system. The circular restricted problem is 
only a first approximation of the problem. The natural problem is the elliptic restricted 


problem. To us it appears that the solution will be of far-reaching consequences 
applicable to solar as well as to stellar system. 


2. COORDINATES OF THE TRIANGULAR LIBRATION PoINTs 
Let us refer our coordinates to Nechvil’s coordinate system (C,& y €) (Duboshin’). 


We shall choose the sum of the two finite masses for the unit mass, the unit of time so 


that the constant of gravitation k2 = 1 and for the unit length, the parameter P of the 


elliptic orbit. Then the equations of motion may be written as 


Boe OH dn 2 OH dG aH 


Si ee 1 
dy Cpe ’ dy ODP , dy Opt | 
apt OH dpe dH dps af () 
‘ BES) as TT Gy eae ae ge 


where 
ic Ay 2 af: 2 2 
= > \ Pe + Pa + Pe) + Pen — pak 


e cos v WwW 
a Rina MIS Dey) 2 R) ok eee ee 
2 (1 + e cos yv) Cr hg 1 + ecosv 


W = (1 — pri + B pro 


a= VEL Ot ae, cee a) Cm 
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#, 1 — » = the masses of the finite bodies (0 < » < }), 


a,p = the reduction factors due to the radiation pressure exerted by the two 
finite bodies, 


(— », 0, 0) = the coordinates of the larger finite body, 
(1 — pu, 0,0) = the coordinates of the smaller finite body. 
For the coordinates of the libration points, we consider the equations 


oH oH oH oH oH 6H 


Tor eg ee Ae 





whence we find that r; = «1/3, ro = B1'3, is a solution of the equations. If L4 and 
Ls be the corresponding libration points, then clearly they are seen to form triangels 


with the finite masses for the other vertices. For convenience, if we puta = SS and 


3 
and B= 5, and (E11, net, Crt, Pez,, Prz,, PS ,,) (i = 4, 5) be taken for the coordinates 
of Li (i = 4, 5), then 


el 
erg Pa Pape em 81 8a.9/P 
s+ 1 — 38, 
eet = 18a V5, Pay = Par Hy bacdeked 


SOS Sa So Ss ee LS ie 0 
2 
2 2 2 92 
Breet — (+9) - 1) [43i+ 3. 


Different from the classical case® here the triangles will be ordinary triangles whose 
sides will be of lengths 51, 52 and 1. 


3. CHARACTERISTIC ROOTS AND THE First ORDER STABILITY 


We shall restrict our study to the planar case alone. Since Ls is symmetrical to 
L4, the nature of motion near Ls will be the same as near L4, so we shall consider the 


motion near Z4 alone. Taking (q¢, pt) (i = 1, 2) for the variations in the coordinates 


of 4, the variational equations may be written as 





peer a die) 088. (19) ...3) 
dy apt’ = dv éqt 
where 
io 2) 


H= > Hk, Ho = constant, H, = 0 
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oy ae es ) ba e cos V ( 2 2 
2 (oi + 3 + pi ge — 41P2 + FOF e cos ¥) ata) 


Eee te. he a: -2 | 
2(1 + e cos »v) 7 {a = 0) 88; tae } 


— _ 3Vvbn |g = ~4 
2(1 + ecos v) [l= 0) 8 82 8, +0 08133"] 


ates [ofa-waetns}—1 | (4) 


] 
1 + ecos Vv 


a 


8 
Cat fw eee = mit 
+ wn (Sn? — 1285) 8,4} + Fah gn Vb {C1 32 88 (582-435) 
7 “ 3 5 
+ 4 513,° (Sy? — 435) } +4n | E87 (1—p) (5b 85 — 1) 


2 2 l Z 
+8, »(5b8) — 1) } + 5H V8 { (1H) 828,765 Sb — 3) 


+ 18,1658; 6 — ah]. (5) 
seniea Ne? lot -6 4 35 
TE aera Pers te ote 


-6 4 15 2 35 
+ ost (ot ths 8 a) 


5 
+3 % % Vb xe - Hw) 808,58 (3-7 Ba," ) 


+ #818;° n( 3-2 923,” ) 


5 3 . 
+ Fah vo{ (1 — 4) £523;° (3—7826 ) 
+ pnd, 353 (3-785 )} 


34 858 ws 
+ 44810 -0ae(- 1+ 58hb + > ea 


35 


ee ghistett 35° ( - EP 
4 Wa ey Bas on b+ 58) +> yn 85° 
pee . 
35 Re an gee qi " 
92 823," b )} - F | 0-0) 3;"( 3-308 32 + 35 2 b2 ) 


(equation continued on p. 407) 
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+p 33°( 3 — 30b 8 + 35 54 b2 )t ..(6) 


Hk = the sum of the terms of the kth degree homogeneous in the variables 41, q2, P1 
and pe 


—E= 58+ 1-8) and y= 5) —1- 8). 


To find the characteristic roots we shall follow Bennet’s method! where he 
expands A, Ho as well as the general solution of the first order variational equation in 
powers of e and also uses 


(1 + ecos v)“1 = | — ecosv + e2 cos? v — e3 cos? v + 


+ (—1)’ e’ cos’ v + ..... eULs 
If the general solution be put as 


y = yO + ye + yl) e 4 ........ ...(8) 


then equating the coefficients of different powers of e on the two sides of the varia- 
tional equation, we get differential equations for y(0), y@), y@) and so on and they can 
be solved successively. Now putting these solutions in the differential equations for 
y(0), y@), »@) and so on, the coefficients A), AG), A@) are calculated. which may be 
given as 


AQ) = 4 [F 4 4{1 — 36n (1 — ») bY”) . (9) 


a{AO)}2 +8 6 
AQ) =0 and A) = TAM }E + 8 A() 


where «, B, y and 5 are constants depending on p. 
In order that the roots (9) may be purely imaginary, it will be necessary that 
pl —-wb<1<S+4b < 0.0285954 (10) 


p = 0.0285954..... corresponds to the resonance case with equal frequencies for b = 1 
and e = 0. This limit coincides with the limit of stability found by Lanzano’. It 
may be noted for distinct 01 and 52 we shall have different ranges of stability given by 
the different values of » for e = 0 shown by us in our earlier paper®, 


Lukyanov® uses Bennet’s method as well as he adopts another method by 
expanding 1/(1 + e cos y) in Fourier series given as 


er csh pe ee (2 « ce'¥ +2 @ 008 204 tee ) 


Vi—e 
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where 


ef = : cos rv d v/(1 + e cos v) 
7 
0 


and then proceeding similar to Bennet he calculates the coefficients A(), X@), A(2) and 


so on, of the charasteristic exponent A. By calculating the characteristic exponents 
by the two methods, he finally compares the ranges of the linear stability, obtained by 
the two methods and he shows that the ranges coincide with those found by Danby?. 


If A) be the characteristic root corresponding to the case « = 0, then 


~ ! 3 1 3 y 
— _ = — e__ — 4— EE 
A +[ (4 mpeg se et Vl—e 


_ 36 "(1 — ub ae 
3 — e 


Ait) 
whence we find that for purely imaginary value of )() 


e < 0.66144...... and (4 — 3/1 — e2)2 > 36 un (1 — p») (1 — e?2). 
yet hi) 
It gives 


we <4 — Vt} — (25— 16 e2 — 24/1 — e2) [36 5. ...(13) 


The equality corresponds to resonance cases with equal frequencies for different values 
of b. We shall not study these cases in further details. The values of » denoted by 
» (e) basing on the range of vales of b varying from 0 to 1 are given in Table I. 


4. NORMALIZATION OF THE HAMILTONIAN FUNCTION Ho 


Taking into view our subsequent studies we shall need the normalization of Ho 
given by (4). For convenience we shall write Ho in the form 


Hz = HY) + AY (14) 


where 
(0) l 2 2 l 9 
H, = z( Pi + P| + Pi 92 — % po + > % (1 — A) 
] 1 9 


HY = (e cos v/2 (1 + e cos v)) ( A gq: +Bargqe+C q, ) ---(14b) 
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TABLE [ 
The values of » (e) 
b e ule) b e u(e) b e ue) 5b e u(e) 


ee tee HE) 2 Pe, Me). 2b ee we) O: 
0.00 00— 


0.7 imaginary 


243 0,0659684 6 00  0,0486645 8 .5 —_0,00753561 
0,05 00—07 - 3 4 0,0429195 1 0.046633 .6 — 0,00139082 
1 00 3 6 0.203579 .2 0,0407785 .7 0,00071481 
ay * = a 0,0037175 3 0,0318208 9 00 —0,0318805 
Sage: 0.376283 3.7 0,0019084 4 0,0209788 1 0,0305737 
ee 0.244755 4 00 0,0750817 5 0,0100732 2 0,026795 
1 4 0.143956 4 1 0,0007185 6 0,00185529 3 0,00209989 
A aes 0,0639155 4 .2 0,0625911 .7 0,00095330 4 0,0138852 
5 Me 0.011237 4 3 0,0485716 .7 00  0,0413961 5 0,00669264 
Ae vscgh 0.005747 4 0,0318206 1 0.03968 .6  0,00123609 
2 00 0 1667 5 0,0151883 2 0,03473 7 0,000635 
got 0 158496 6 0,0278553 3 0.02714 1 00 0,0285954 
2 0.135784 7 0,00143063 .4_—(0,017918 1 0,0274273 
Ka 01030426 .5 00 0,0590414 5  0,00862 2 0,0240476 
Ce 0,0659679 0,0565477 6  0,00159 3 0,0188399 
ee 0,0308681 5 0 0493768 .7 0,000817 4 0,0124789 
2 6 0,0058676 3 0,0384482 .8 00 —0.0360196 5 0,0060192 
a 4 0,0028654 4 0,0252858 1 0,0345365 6 — 0,0011123 
3 00 0.1032539 5 0,0121128 .2 0,0302519 7 0,000571 
3 J 0,0986483 6 0,0022272 3 0,0236663 
1 2 0.0853556 —.7 0,00114428 4  0,0156489 


3 fa—pest+ uns? 


} 

A= 

| 
p= Svb{eea-w+iaeft (15) 

| 

2 
can fa-—ogt+ nah. | 
For normalisation we shall introduce the transformation 

-.-(16) 


(41, q2, Pi, p2) = (4; ’ I : Py» Ps ) N : 
where N is the same as used in Manju and Choudhry?®. 
Here 


2 (0) )? 
oF = — {a0} and os re - {ani} 
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we shall cal] ; and we as the frequencies and they are given by the equation 
wo4— w2 + 9u (1 — p) b= 0. <as(17) 


The transformation (16) reduces the Hamiltonian (14) to-the form 
1 ] 12 2 12 1 12 2 12 
Hi = } (ep? + of a?) - 5 (a + 0 a2") 


e cos V g Spe 
+ T+ ecosy >, [Pod aie oe OSM 


utt=2 


where v = (v1, v2), » = (1, v2) and v1, v2, wy and pe so vary that v + p = vy + ve 
+ 1 + »1 = 2 and for simplicity we shall not write the ranges again and the coeffi- 
cients a‘, are given as follows : 


‘ 


bis 2 j 1 ¢ 2 
25000 = 5 (A+ Bat Cc), Ayo00 = 74 ( A+ Ba+Ce} ) 


' Pai by Bed ee l 2 12 
M00 = > © 41 91s A092 = 7 ©% 5, 


‘ ] 

21499 = 41 42 {4+ 7 B(cy + co) + Care 

25010 = a, by (§ B+ Cy), 24001 = — a, a2 be (4 B+ Co}) 
A119 = 91 92 51 (4 B+ Cee), ao,, = — a1 43d, bo C 


I%01 = — a; be (4 B+ Ces). 


Next we shall introduce the following transformations : 


~ 
‘ 


l ‘ a 
 mley real Wt (Moa e Ao et AC Le) 2) 


) 


S (0 = 
which transform Hi to Hy” given as 


A, = harp, +a) — too (p +43 ) ...(20) 
WY Ad 
and H, to H, given as 


C8 Od Calls fad pea 


(i)_ se Cos Reet Sek en 
as 1 + ecos vy 3 Wed) Ga Pra ike -- (21) 


uth? 
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where ayy can be easily calculated. For the convenience of subsequent calculations 
we shall introduce the complex conjugate canonic variables given as 


~~ 


q, =Pitigqin Py =pi—ig (i =1, 2). ...(22) 


Consequently the Hamiltonian will be reduced to the form H; = 2i He 


where 


Hy = iotq; Py — ies PS 


.  e cosy ee a re 
a. 2i l oo ecos v Ss Gy gy BY Py Py tea) 
vtp=2 
The coefficients in H} aresuchthat a’, = aj 
v 
where the bar sign denotes the complex conjugate quantity. Other coefficients are 
given as follows : 


~ 


Bso99 = 4 (— 42000 + ao020 — i 41010) 


~~ 


Ayo000 = *(— 40200 + Go002 — 7 Go101) 


~ ~~ ~ 


Beinn = }(- ai100 + 40011 — 7 41001 — i 40110) 


tai = }(— ai100 + 0011 — 7 41001 — 7 ao110) 
Qso10 = + (42000 + 40020), 4q49, = 4 (40200 + ooo2). ...(24) 


Next we shall find the transformation 


( 9; Ps ) os (:eF". ps" ) (25) 
reducing the Hamiltonian (23) to the normal form in complex conjugate variables 
given as 


Hy Cre Pon ) = id g** py* + i A293" PS”: ...(26) 
Let this transformation be given by the generating function 


~ ‘ F * *e 
q3 prt + P59, S {ai a5 <P » Po v) 
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where 


"Uy kek 
; i 2 1 Fg 
S= Syu q q, Py Py 


vt+Ph=2 


and sy, are to be chosen 27-periodic functions of v. The relation between the vari- 
ables g;, p; and gj", p?* are given as 





as a es 
FMCG fos merentt ys Se 
ap** éps 


whence we the identity 





és Os Os 
H3* ( 95 + Pv )— Hs ( af, of +e) = 
ép;* op; 


On expanding and equating the coefficients of equal powers on the two sides, we shall 
get 





























2 oH** 
; as 2 és \2 62H. 
Hy" ( 45.23.) + ke i ee 
j=1 OP; 0q; op, qi 
e2 H** o2 H** 
0 é 2 2 2 
sO rc ecole (I ) ] 
op, Org Oa, Ogg: | aps® “ie age 
2 @ H* 
3 2 
=; ( ejapten) > 
- g-4 OQ; Ape” 
02 HY’ 
as \2 2 
ape kg) Neer t sei 2 
“ ** « a“ 
ns eq; 693 
62 H* eg Ps 








é 2 
<a | eee 
opt* apx* 2% opx*? 


sion. chee sale aden 


2 


yth=2 
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whence we get 


413 
. « . a“ . x te 
Pargie, 4+ 3 Aad,-Ps) ti > (ur Ar + po As) Sve Gy” 


u+ P=? 


° « ,** : w 
—iayg; p, + iw. 9, po 
.  ecosy 
—2 


"y "y eK tig 
a’ q Ei q 2 Pp gt p 2 
1+ ecosy yr AY 2 1 2 
V+p—2 


chy er eh) **P 
(v1 1 — v2 @2) Syn 9, Ye 


2 
Py Py 


V+p-2 


"y ” 
ep dsyp See 


xe 
« 1 2 
eet. ae 
vtPh=—2 


Restricting only upto the second order terms in e, we may write (28) as 


irq, PL + ide gape! SS (v1 At + He Ag) [e s@) + e? s)) 


. . . 2 
— 101 q; Tee + i 2 gPar hae F | € cos y— 5 + €08 2») | 


Se 
“ 
4 > Ayn 4, I, Py Po 


— i & (vy &1 — ve 2) (e sQ) + e2 s(2)) 








ds (1) ds (2) 
a 2 
fen. ce ...(29) 
where 
Syn = e yy st) +e 3?) ; 
Let us write 
} 
et ea + et ay + 
and r ...(30) 
eee ere pe 
J 
and equate t 


he coefficients of the equal powers in e, then we shall have 


3 0 
a” hy = we, EAD + pee +i Sy we hye + pe ai ) ) s@ 


2 
ds() 
mot, Of COS V Ay, Pos | 3 (v1 @1 — V2 w9) = ie ...(31)) 


1 G1, 
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whence we shall have the following relations. 





As? } 
1010 aa) (0) (1) Peay 
gti im GAP HEAL Siorq — 24 COSY Ajoiq — #1 Si010 | 
| 
(1) | 
ds 
0101 ; F : 3 1) > 
Sa * lala i Ay +i Me Aes — 21 COSY Aj19, — I We hol | 
1 
ds* ae | 
= + i [(1 — p1) 1 — (ve — pe) 2] sD = — Zicos va‘, } 
£32) 
On integration, we shall have 
jak = eye v — 2 SiN V A610 
iat = — 2isinvayi9, + ia ! 
(2) _ 2é a’vm [sin v + i (v1 — 41) w1 — (v2 — M2) 2} cos v) (33) 
Svps [(v1 — 1) ©1 — (v2 — pe) we]? — 1 oa 
(1) (1) 


(1) (1) 
AV =A, 


By virtue of periodicity of s},,,. ands), it follows that = 0 and the re- 


lations (33) completely determine S as a complex-valued function restricted to the first 
order terms in e alone. 
Now let us find a real-valued transformation (qj, pj) > (q;, P; ) reducing 
haw ae 0 salt} 
the Hamiltonian Hz = Hi + H, tothe normal form given as 


sie 4A ( e” *° } «? e” 
— 1 q, + P,; )+ 9 Ag (9, + Ps ) -..(34) 


~~ 


Let this transformation be given by means of the generating function q P; + qe 


* werle : : 
Po + K(qis DP; v), where K is restricted to the order of e alone. From the trans- 
formation formula 








3 he OK 
q; I= Pp, - + a? 
ae J 
a p, oq, 


we may obtain by implicit function theorem and using the fact that all the variables 
are small 


~ ‘ @K ~ . OK 





: (35) 
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where K* will be obviously of the order of e. 
By the formula (26) correct to the order of e we have 
: os** ; as** 
eee ee 
é P; 0 q; 


where S** = S(1) (qg**, p**, v). Further taking into account the relation between 
the complex canonic variables with the real ones 


...(36) 


q, = Pit+igqi Pp; = pi -i Qi 


a * 


eee, Ee hy =P, 19; 


and denoting the function S() (p* + i q*, p* — ig*, v) by W (q*, p*,v), we shall 
obtain by the formula (36), 





es 1 OW - le <0: 
a t,t 5 ..(37) 
é P; 7] q; 
Comparing (35) and (37), we obtain K* = iw... ...(38) 


and the function K= Ky, q; J .° pr Ps a 


the formulae (4b), (19) and (33) its coefficients may be found which are given as follows. 


will be realvalued. With the help of 











K2000 = =-( = feed a ee 23 ee ) | 

ko200 = ae ( 3 eee - Chee BS ae ) 

koo20 = sr ne uP Gan 3s er ) 

koooz = + fe a5 Poi ar he ) 

k3100 = s( a + Scot $5 ee or et ) 

kino = she i: ae 

kioo1 = 2 are 34 oni st ae i Spon 

koi10 = 4 ( eee Ai ‘Sl + cee z spon 

kool = leek = ake | 

koou1 = a ( Hee + pe ne Hits e sn ). ..(93) 
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Thus the transfarmation of the Hamiltonian H2 to the normal form given by (34) cor- 
rect to the first order of eccentricity has been found. It is obtained through the trans- 
formations (16), (19) and (35) and the coefficients of the generating function K are 


given by (39). 


5. RESONANCE CASES 


Taking into view the complete study of the stability of the triangular libration 
points we shall need an investigation if resonances are present. Since we aim to apply 
KAM-theorem for stability which needs mainly the third and the fourth order terms, 
we shall only examine it the resonances of the third and the fourth order exist. 

If a curve giving resonances be plotted for different values of e, then we shall 


need the value of A; and Az at least correct toO (e), since an = pe = 0. The qua- 
ntities A and ye are found by the periodicity conditions of the functions Phil and 


Ke For it let us take up the expansion (27) and here equating the coefficients 


with e?, we shall get 


s 





(2) 
x ae = — 2icosv| 4a (1) +a‘ (1) 
pe 0020 *2000 1010 1010 
” (1) « (1) - ; : 
+ 4599; $3001 + 20110 i100 ) + 24 £08" Y Qor9 
~~ (2) 
+ iA, 
(2) 
4S o101 by (1) (1) (1) 
F tale Mai aS (4 Aoo2 Soz00 + 4110 51001 + 70101 50101 


“ (1) . - te 
+ Geass wre | + 2i cos? va5,9, + as ) 


Substituting into the right-hand side expressions of these equations the value of the 
; ees : 
functions s(t given by (33) and choosing he and Ye such that the constant terms 


on the right-hand side are equal to zero (the condition of periodicity of ae and 


s®) ) we shall obtai ipulati 
0101) We Shall obtain after some manipulations using the formulae (14b), (18) and 
the equation 


o4 — w2 + 9un(1 — pbb = 0 


the following expressions for ve! and aye and obtained. 
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@1 2 (6 wt — 7) 


os 
2 2 ; 
4 (4w, oe 1) (2; —1) 
2 2 
@ 7 @2 w) (6 ig PY 7) 
2 ...(40) 


et 
4(40, —1)Q 0 ~1) 


These values coincide in forms with Markeev’s!9 values. 


Let the value of yw giving the resonance ky Ai + ke Ap = N for small e and 
correct to O (e?) be given as » = (9) + e2 u(2) where #() is the value of » when 
e = 0 and x (?) is the contribution to the value of » when e ~ O and it is considered 
correct to O(e?). Letting Ar = Aq (u) + e2 w(2)) and Az = Ao (u() + e2 u(2)), we 
have on expansion by Taylor’s theorem that 

(0) (2) dr 
a emtX 2 2) 
Tat Feta, be ne (3 ), 
and similarly, 
(0) (2) 
Ao =A, tera, + e2 nO) (3) 
2 a du Jo 
(2) (2) ; 
where A, and A, are given by (40) and the suffix (0) denotes that the value is to be 
taken when » = (9). Putting these values in ky Ay + ko A2 = N and equating the 
coefficient of e2 to zero, we shall have 


(2) (2) 








ri(2) = : 
k dog = k dws 
2 dp 1 dp 


The value of (2) is calculated on putting » = »() on the right-hand side. The values 
of »() and p(2) for the different resonances of the third order are given in Tables II 
and III. In Table IV we have denoted by » (e”) the values of » giving the limit for the 
linear stability for different values of b and e under the resonance case 3A2 = — I. 
The values have been calculated on the formula »p = »(°) + e? »@), e has been taken 
to vary frome = 0.0 toe = 0.7. We have here compared the values of u (e) and 
uw (e2) for the resonance under reference and we have examined if the resonance exists 
for a given value of b and e and if the study of stability will be necessary in the parti- 


cular resonance case when higher order variational terms are taken into consideration. 


The conclusions are listed in the succeeding section. 
We shall examine the following resonances of the third order if they exist under 
our range of linear stability given by Table 3.1 : 
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TABLE II 
The value of (°) 


TE 











atic | ees ie Deas Dry +r. = 1 Ay — 2%, =2 3, = —2 
b pl?) 2°) (°) »(°) (9) 
Ai 0.1255 0.23126 imaginary imaginary imaginary 
2 0.0583 0.09861 0.15071 0.15071 0.1641 
3 0.03803 0.06326 0.09421 0.09421 0.1018 
4 0.02823 0.04662 0.06872 0.06872 0.07407 
5 0.02245 0.03692 0.05413 0.0513 0.05826 
6  0,01864 0.03056 0.04466 0.04466 0.04303 
7 0.01593 0.02608 0.03802 0.03802 0.04086 
8  =0,01391 0.02274 0.033095 0.033095 0.03556 
9 0.01235 0.02016 0.029303 0.029303 0.03147 
1 0.0111 0.01811 0 026291 0.026291 0.02823 
.75 0.014853 0 24294 0.0353854 0.0353854 0.038026 
w(%) pl) = 4 — Nee SR = ne yO) = rf a Wiel S sep = oo fa — 
ae ee ra gay MO cr 2 Bp 
TaBLe III 
The values of w\*) 
2. = — 1 Ay + 2A. =O 2h, +>, = 1 Ay — 2g = 2 3, = —2 
b 
a —0.83513 —3.8097591 imaginary imaginary imaginary 
2 —0.35404 —1.2733528 0.51445 1.5010664 0.568121 
3 —0,22567 — 0.780191 0.295214 ().8613774 0.3198794 
4 —0.165736 —0.56367 0.2083243 06078504 0,2242903 
Me — 0.130984 —0.441488 0,1612059 0.470368 0.1730103 
6  —0,10829 —0.3629225 0.1315443 0.3838212 0.140912 
or —0,09290 — 0.308136 0.111132 0.3242611 0.1188955 
8, —0.080427 —0,267732 0.9962145 0,2807356 0,1028464 
9,  —0,071262 — 0.236704 0.084835 0.2475324 0.09621 
1 —0,063972 —0,212128 0.075866 02213626 0.0809987 
we) = en u(2) xo 220764256 _ 
(1 — 2p) = 6 A — 2p) 


eee 
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(i) 3A2 = —1, (ii) Ar + 2A2 = 0, (iii) 2A2 + Ag = 1 
(iv) Ai — 2A2 = 2, (v) 3A2 = —2. 


’ 


The corresponding values are given by Tables 1V—VIII. 


Coming to resonances of the fourth order we may proceed similarly and it may 
be seen that the study of the stability will be needed for the following resonance. 


(i) 44g = -1, (ii) Ay + 3A2 = 0, (iii) Ar — 3A2 = 2, 
(iv) 2 (Ay + Ag) = l, (v) 3A1 + As = 2, (vi) 3Ay — Ag = - 
(vii) Ar + 3A2 = —1, (viii) 4A, = 3. 
TABLE LV 
The values for  (e*) 


(3,2 = — 1) 
b e = 00 e= 0) e= 072 ea= 0:3 e= 04 eo U:5 e = 0.6 e= 0,7 


— MG <iscseae ee = 


Ot.) *O.12355 0.117 0.092 0.05 negative —ve —ve 














0.2 0.058 0.055 0.044 0.026 0.002 —ve —ve as 
0.3 0,038 0.036 0.029 0.018 0,002 —ve —ve —ve 
0.4 0,028 0,027 0.022 0.013 0.002 —ve —ve —ve 
0.5 0,022 0,021 0.017 0.011 0,001 —ve —ve —ve 
0.6 0.019 0.018 0.014 0.009 0.001 —ve —ve —ve 
0.7 0,016 0.015 0,012 3,008 0,001 —ve —ve —ve 
0.8 0.014 0.013 0.011 0,007 0,001 —ve —ve —ve 
0.9 0.012 0.012 0.009 0.006 0,001 —ve —ve —-ve 
1.0 0,012 0.010 0,009 0,005 0,001 —ve —ve —ve 
—ve = negative 
TABLE V 
The values for » (e*) 
(A, + 2dr. = 0) 

b e=0,0 e= 0,1 c = 0,2 e=03 €=04 e¢=05 e=06 e=08 
oa O23 0,193 0,079 —ve —ve —ve —ve —ve 
0.2 0,099 0.0859 0.048 —ve —ve —ve —ve —ve 
0.3 0,063 0,055 0,032 —ve —ve —ve —ve —ve 
0.4 0,047 0,041 0.024 —ve —ve —ve —ve —ve 
0.5 0.037 0,033 0,019 —ve —ve —ve —ve —ve 
0.6 0,031 0,027 0.016 —ve —ve —ve —ve —ve 
0.7 0.026 0,023 0,014 —ve —ve —ve —ve —ve 
0.8 0,023 0,020 0.012 —ve —ve —ve —ve —ve 
0.9 0.020 0.018 0,017 —ve —ve —ve —ve —ve 
10 0.018 0,016 0,001 —ve —ve —ve —ve —ve 


SO 
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As shown by Moser2, the motion under the resonances A, — 3A2 = 2 and 3A, —A2=3 
cannot lead to instability and the cases Ay + 3A2 = O and Ai + 3A2 = — : reduce 
to the case Ay + 3A. = m (an integer). Thus it remains to study the stability only 
for the following five cases : 


(i) 442 =m, (ii) Ar + 3A2 = m, (iti) 2 (Ar + Ao) = m, 
(iv) 3A, + Ag =m, (v) 4A. =m, 


where mm is an integer. 








TABLE VI 
The values for v- (e*) 
(2A1 + 4s = 1) 

b e=00 e=0,1 = 0,2 = 0,3 e= 0.4 e=05 e=0.6 e = 0.7 
0.1 Imag. Imag. Imog. Imag. Imag. Imag. Imag. Imag. 
0.2 0,151 0.156 0.017 0,197 0.233 0,279 0,339 0,403 
0.3 0.094 0.097 0.106 0.121 0.141 0,168 0,200 0,239 
0.4 0.069 0.071 0.077 0,087 0.102 0,121 0.144 0.171 
0.5 0,054 0.056 0.061 0.686 0.080 0,094 0,112 0,133 
0.6 0,045 0 046 0.050 0,066 0,066 0.078 0,092 0,109 
0.7 0,038 0,039 0.042 0,048 0,056 ),066 0,078 0.092 
0.8 0,033 0.034 0,037 0,042 0,048 0.057 0.068 0,081 
0.9 0,029 0.030 0,033 0.037 0.043 0.051 0.060 0,070 
10 0,026 0.027 0,029 0.033 0,035 0.045 0.054 0,063 

TABLE VII 


The values for » (e") 
(A; — 2A, = 2) 


b e=00 e=0.1 e=02 e=03 e=04 e=05 e = 0.6 e= 0,7 
0.1 Imag. Imag. Imag. Imag. Imag. ; Imag. Imag. Imag. 
0.2 0,151 0.166 0.211 0.286 0,391 0.526 0,691 0.886 
0.3 0,094 0,103 0,129 0.172 0.232 0,340 0.404 0.516 
0.4 0.069 0.0748 0,093 0,123 0,166 0.221 0,288 0,367 
0.5 0,054 0.0588 0.0729 0,0965 0.129 0.1717 0.2235 0.2846 
0.6 0,045 0,048 0,060 0,079 0,106 0,141 0,183 0.233 
0,7 = 0,038 0.0413 0.051 0.0672 0,09 0.1191 0.1548 0.1969 
0.8 0,033 0.0368 0,0443 0.0623 0,085 0.1142 0.1498 0,192 
0.9 0,029 0.032 0.0392 0.0515 0 0689 0.0912 0.1184 0.1506 
10 0,026 0,0285 0.035 0,046 0.062 0,082 . 0,106 0.1348 


Imag. = Imaginary 
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TABLE VIII 


The values for u (e*) 





3==—2 
0.1. Imag. Imag. Imag. Imag. Imag. Imag. Imag. Imag. 
0.2 0,104 0.170 0.187 0,215 0.255 0.306 0,369 0.443 
0.3 0,102 0.105 0.115 0.131 0.153 0,182 0,217 0,259 
0.4 0.074 0.077 0.083 0,094 0.110 0.130 0.155 0.184 
0.5 0.058 0.060 0,065 0.074 0,086 0,102 0,121 0.143 
0.6 0,048 0,049 0.054 0,061 0,071 0,083 0,092 0.117 
0.7 0,041 0.042 0.046 0.052 0,060 0,071 0.084 0,099 
0.8 0.036 0,037 0.040 0,045 0.052 0,061 0,072 0.086 
0.9 0,031 0,032 0.035 0,040 0.047 0.056 0,066 0,079 
1.0 0,028 0.029 0.331 0.036 0.041 0.048 0.057 0,068 


OS ee eee 


Imag. = Imaginary 


6. CONCLUSIONS 


In section 3 we found out the characteristic exponents correct to O (e2) and we 
have calculated the values of » giving the range of linear stability for different values 
of e. These values denoted by » (e) have been given by Table I. In section 4 we have 
normalized the second order terms by Birkhoff’s transformations and also we have 
calculated the various terms as resulting after transformations. In section 5 we have 
found the values of » corresponding to different types of resonances. e has been taken 
to vary from 0.0 to 0.7 and 5 from 0.1 to 1.0. From Table IV we find that the 
resonances of the type of 3A2= —1 exists fore < 0.5 and in all such cases the motion 
is stable in the linear sense since the corresponding values of » (e) is greater than 
p (e?). So for the investigation of non-linear stability this type of resonance case has 
to be taken into account. The negative values of » (e2) indicate that the resonance 
case will not exist for our range of values of » (e) under consideration. 


Table V shows the resonance case Ai + 22 = O exists for b = 0.1 only for 
values of e such that 0.1 << e < 0.3, and for b > 0.2 it exists fore < 0.2. Since the 
values of » (e2) corresponding to the resonance are all less than the value + (e) cor- 
responding to the range of linear stability, for non-linear stability the study of the 


resonance case will be necessary. 


Table VI shows that the resonance case 2A1 + Ag = 1 will not exist for b = 0.1 
- 0.2, the resonance exists for all e from 0.0 to 0.7, but the the linear 


and for b2 
> 0.2 the motion will be unstable 


stability is possible only fore < 0.2, and for e 
since even the linear stability does not hold. 
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Table VII shows that for b = 0.1 the resonance Ai — 242 = 2 does not exist 
but for values of b it does exist. The values of » (e2) are greater than the correspond- 
ing values » (e) for the linear stability. [t shows that the motion can be stable only 
for the values of e such that e < 0.1, and for e > 0.1 the motion will be unstable. 


Table VIII shows that in the resonance case 3A2 = —2, the motion can be 
stable only for e < 0.1, and fore > 0.1 the motion will be unstable since p (e) 
< p (e?) for all b and all e. 


Similarly resonances of the fourth order will exist for some values of b and e 
and for non-linear stability these cases have to be taken into account. We have not 
calculated the corresponding values of » (e”), but it is guessed that the resonances of 
the fourth type will exist, similar to the classical case as examined by Markeev!9. 
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A paper entitled “Steady Plane MHD Flows with Constant Speed Along Each 
Streamline” authored by M. A. Sattar was published in Indian J. pure appl. Math. 
18 (1987), 548-56. Exactly the same work authored by M. A. Sattar and O. P. Chandna 
appeared in J. Math, Phy. Sci. 22 (1988), 321-33. Although Dr Sattar did this work 

jointly with me, I was totally unaware of its submission to these journals by him. The 
mathematical analysis of this work has two cases when 
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is nowhere zero in the region of flow and when J is everywhere zero in the region of 
flow. The summing up of the case when J = 0 everywhere needs correcting. This sum- 
ming up should read: The steady plane flow of a viscous incompressible fluid of in- 
finite electrical conductivity having constant velocity magnitude along each streamline 
has only constant solutions given by (3.31) when J = 0 and the streamlines for such a 
flow are parallel straight lines. 


The above correction implies that the non-constant solutions given in (3.35) are 


not possible for the subcase taken. This subcase has 


u2 + ¢2(u) = u2 4+ v2 = constant = c-say 


ux + $' (u) Wy = 9 ...(3.27) 

fu2 + ud (u) $! (u)} uz + {ud (u) + ¢ Ww) ' (u)} uy = 0 ...(3.28) 
(2u  (u) + 62 (u) ou) — wo! (Wh Ue 

...(3.29) 


+ {$2 (u) — v2 — 2ud (u) ¢' (u)} uy = 0 


u 
where c is an arbitrary constant and = x 0, 
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Since 2u+2¢ (u) 4’ (u)=0 is a consequence of the assumption u2 + ¢2(u) = con- 
stant taken as Case III in paper, it follows that equation (3.28) is identically satisfied. 





a 
> from equations (3.27) and (3.29) and requiring a + 0, we find 
that ¢(u) must satisfy 
¢2 (u) 6'2 (u) — w2 62 (u) — $2 (u)+ u2 + 4ud (u) ¢' (uv) = O. 


Eliminating 


Employing ¢(u) = + Wc — wu? in this equation to be satisfied by ¢ (u), we conclude 
that 


c= 0 
and, therefore, 
u2 + $2 (u) = u2 + v2 = 0. 


This equation implies that subcase III], for the case J = 0, given in the paper is not 
physically possible. 
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